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Abstract 

The beta decay of molecular tritium currently provides the highest sensitivity in laboratory-based 
neutrino mass measurements. The upcoming Karlsruhe Tritium Neutrino (KATRIN) experiment 
will improve the sensitivity to 0.2 eV, making a percent-level quantitative understanding of molec¬ 
ular effects essential. The modern theoretical calculations available for neutrino-mass experiments 
agree with spectroscopic data. Moreover, when neutrino-mass experiments performed in the 1980s 
with gaseous tritium are re-evaluated using these modern calculations, the extracted neutrino mass- 
squared values are consistent with zero instead of being significantly negative. On the other hand, 
the calculated molecular final-state branching ratios are in tension with dissociation experiments 
performed in the 1950s. We re-examine the theory of the final-state spectrum of molecular tritium 
decay and its effect on the determination of the neutrino mass, with an emphasis on the role of the 
vibrational- and rotational-state distribution in the ground electronic state. General features can 
be reproduced quantitatively from considerations of kinematics and zero-point motion. We sum¬ 
marize the status of validation efforts and suggest means for resolving the apparent discrepancy in 
dissociation rates. 
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I. INTRODUCTION 

The fact that neutrinos have mass [1, 2] is the first definitive disagreement with the 
minimal Standard Model of particle physics. As new extensions to the model are developed, 
a determination of the absolute neutrino mass scale will be essential [3]. In addition, this 
mass scale influences the large-scale structure of the universe and is an important ingredient 
in cosmological models [4, 5]. Observables related to the neutrino mass are accessible through 
cosmological studies, neutrinoless double beta decay, and supernova neutrino observations. 
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However, the most direct approach to the neutrino mass, with minimal model dependence, is 
by detailed measurement of the shape of the nuclear beta-decay spectrum near the endpoint. 

Tritium (T) undergoes an allowed nuclear beta decay, transforming to 3 He with the 
emission of a beta electron and electron antineutrino. The low Q-value of 18.6 keV means 
that the modification of the spectral shape by the neutrino mass is relatively large. In 
addition the half-life of 12.3 years allows sources with high specific activity to be constructed. 

The well-known form of the tritium beta spectrum is illustrated schematically in Fig. 1 
for massless neutrinos and for 1-eV neutrinos. It is the task of the experimentalist to 




Kinetic Energy - Q-value (eV) 


FIG. 1. Tritium beta spectrum with 3 active neutrinos with masses m U i ~ 1 eV for the case of no 
daughter excitation. The left panel shows the full spectrum. The right panel shows the last 5 eV 
before the endpoint, with the dotted curve indicating the spectral shape for = 0. 

measure the spectral shape and thereby determine the neutrino mass. Only a fraction of 
order 10“ 13 of the decays populate the last 1 eV of the beta spectrum. Uncertainty on the 
Q-value and practical experimental challenges preclude fixing the endpoint energy during 
data analysis and it is therefore treated as a fitted ‘nuisance parameter.’ Furthermore, the 
spectral distortion due to the neutrino mass is small and distortions of similar size can arise 
from a number of theoretical corrections and from instrumental effects. For a molecular 
tritium source, the largest modifications to the spectrum are caused by excitations of the 
daughter molecule formed in the decay, which must be calculated from theory. One could 
consider using a non-molecular source, such as T + or T, but these are far less practical due 
to space-charge limitations and the high reactivity of atomic hydrogen. 

The ongoing construction of the Karlsruhe Tritium Neutrino experiment (KATRIN) [6], 
the next-generation tritium-based neutrino-mass experiment, has renewed interest in the 
molecular final-state distribution (FSD) populated by T 2 beta decay [7]. With a design 
neutrino-mass sensitivity of 0.2 eV, KATRIN depends critically on a theoretical under¬ 
standing of molecular effects. Accordingly, extremely precise, ab initio calculations of the 
molecular final-state spectrum have been performed [8, 9] in the region of interest for KA¬ 
TRIN, near the endpoint of the beta electron energy spectrum. A direct experimental 
verification of these calculations through a study of the molecular final-state spectrum itself 
is not practical as explained in Sec. VI. Indirect tests can be performed, but have yielded 
mixed results. Although most of the spectrum of the HeH + isotopolog is inaccessible to 
experiment, many predicted spectral features have been observed in emission; HeH + pho¬ 
todissociation measurements are also compatible with theory, although a high-precision test 
has yet to be performed. On the other hand, measurements of the branching ratio of T 2 to 


3 




the bound molecular ion 3 HeT + following beta decay - another observable indirectly related 
to the final-state distribution - show stark disagreement with predictions. 

In this work, we discuss the aspects of the neutrino-mass measurements that have moti¬ 
vated study of the molecular final states excited in To beta decay and summarize the current 
state of theoretical work on the topic. We begin by examining the commonly used theoreti¬ 
cal expression for the spectrum of allowed beta decay and derive a more general expression 
that facilitates a consistent treatment of molecular ‘final-state’ excitations. Focusing on 
the region of the spectrum near the endpoint, we show that the energy spread caused by 
molecular excitations is dominated by the zero-point motion of the parent T 2 molecule, and 
derive a general analytic expression for the variance of the ground-state manifold of states 
that includes not only zero-point vibration, but rotational and translational degrees of free¬ 
dom. The expression can be applied at any selected temperature up to 300 K, to the 3 
isotopologs T 2 , DT, and HT, to any chosen ortho-para admixture, and to address possible 
uncertainty in the rotational-state temperature. The variance of the final state distribution 
is found to be quite sensitive to whether rotational thermal equilibrium has been achieved 
in the source gas. We then examine several indirect experimental approaches for validating 
theoretical calculations of the final-state distribution, and review existing measurements. 
When modern calculations are used to re-evaluate gaseous tritium experiments performed 
in the 1980s, it is found that negative values of ml are eliminated. We suggest desiderata for 
a new experimental investigation of the branching ratio to the ground-state manifold with 
a view to resolving the discrepancies of more than 50 years’ standing. 


II. DIRECT NEUTRINO MASS MEASUREMENTS: EXPERIMENTAL PROGRESS 

A. Historical tritium-based neutrino-mass experiments 

Tritium-based experiments to measure the absolute mass of the neutrino have a long 
history. Robertson and Knapp [10] review early experiments, Otten and Weinhcimer [7] give 
a detailed treatment of more recent experiments, and Drexlin et al. [11] review experiments 
that are currently under construction. 

The issue of atomic and molecular excitations in tritium-based neutrino experiments was 
first raised by Bergkvist in the early 1970s [12]. He was able to set a 55-eV limit [13] and 
noted that an understanding of daughter excitations was required to improve limits further. 
His work motivated the construction of an experiment with a windowless, gaseous T 2 source 
at Los Alamos National Laboratory (LANL) [14, 15]. The use of T 2 is advantageous because 
the molecular final-state calculations are more tractable than for more complex sources, and 
a gaseous source minimizes the effects of scattering on the beta spectrum. The LANL 
experiment yielded an upper limit of m v < 9.3 eV at the 95% confidence level [15] with a 
2-cr excess of events observed in the endpoint region, reported quantitatively as a negative 
central value of ml. An experiment at Lawrence Livermore National Laboratory (LLNL), 
also using a windowless, gaseous T 2 source, yielded a central value in good agreement with 
the LANL result, but with much reduced statistical uncertainties. The excess of events near 
the endpoint then corresponded to 6cr [16]. 

Concurrent experiments in Beijing [17], Tokyo [18] and Zurich [19] used complex tritium 
sources. All of these experiments gave results that were consistent with zero neutrino mass 
but with central values in the unphysical negative-mass-squared region, which is symp¬ 
tomatic of an underestimated theoretical or experimental contribution to the resolution 
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function. Attempts to reduce such influences furthered interest in molecular-tritium experi¬ 
ments, where ab initio molecular calculations were possible, and inspired further theoretical 
work on the molecular final-state distribution in the late 1990s (Sec. IV). 

The Particle Data Group evaluation [20] of the present limit on the neutrino mass, m v < 
2 eV at an unstated confidence level, is derived from the Mainz [21] and Troitsk [22, 23] 
experiments, both of which employed a new type of spectrometer. In a magnetic-adiabatic- 
collimation-with-electrostatic (MAC-E) filter [24], the momenta of beta electrons rotate to 
a mostly longitudinal direction as the electrons pass from a region of large magnetic field 
to a region of magnetic-field minimum. The kinetic energy of the resulting broad electron 
beam is then analyzed with a longitudinal retarding potential. 

The Mainz source consisted of T 2 films quench-condensed onto substrates of highly ori¬ 
ented pyrolytic graphite. Solid-state source effects, such as dewetting effects and local lattice 
relaxation after the decay of a bound tritium atom, required careful attention in the Mainz 
analysis. The final Mainz result was m u < 2.3 eV at 95% confidence [21]. 

The Troitsk experiment, like its predecessors at LANL and LLNL, used a windowless, 
gaseous tritium source. The gas density and source purity were monitored indirectly by a 
mass analyzer at the source and by count-rate measurements at a low retarding-potential 
setting. During later runs an electron gun mounted upstream of the source was used to 
monitor the column density. The initial analysis of the data required the inclusion of a 
step function added to the spectral shape [22], the so-called “Troitsk anomaly.” The final 
Troitsk result, based on a re-analysis of the subset of runs for which electron-gun source- 
column-density calibrations were available, was rn u < 2.05 eV at 95% confidence [23]. No 
step anomaly was required in the re-analysis. 


B. Future prospects for direct neutrino-mass experiments 

As the sensitivity of T 2 -based experiments improves, an accurate understanding of the role 
of molecular final states after beta decay becomes increasingly important. The systematic 
uncertainty associated with final states has been a major motivator in the search for other 
experimental approaches to direct neutrino mass measurement. The common alternative 
approach employs microcalorimeters with sources of rhenium (MARE [25]) or holmium 
(HOLMES [26], ECHo [27, 28], and a LANL experiment [29]). Microcalorimeters suffer 
from pile-up spectral distortions, requiring the construction of a large number (order of 
millions) of functionally identical calorimeters. 

Alternative measurement techniques using tritium sources are also being explored. An 
approach for coincidence detection of the beta electron and the 3 He + ion from a source of 
trapped tritium atoms was proposed [30] but later shown to be infeasible [31, 32], The 
Project 8 collaboration is currently studying the feasibility of measuring beta electron en¬ 
ergies by trapping and measuring their cyclotron radiation frequencies with microwave an¬ 
tennae [33, 34], In its planned use of a T 2 source, Project 8 again requires knowledge of the 
molecular final states of the source, although the collaboration is also studying the possibil¬ 
ity of building an atomic T source by magnetically trapping single atoms as well as emitted 
electrons. Substantial research and development are required before a full experimental 
design can be developed. 

Molecular-tritium beta decay remains the major focus of experimental work on the direct 
measurement of the neutrino mass. Scheduled to begin taking data in 2016, the KATRIN 
experiment will be the most sensitive neutrino mass experiment to date with a design- 
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sensitivity of 0.2 eV at the 90% confidence level [6]. To achieve this level of sensitivity, 
the total systematic uncertainty must be controlled to within a budget of approximately 
v sys t{ml) ~ 17 x 10” 3 eV 2 . 

The molecular final-state distribution populated by T 2 decay represents one of the larger 
potential sources of systematic error in KATRIN. A 1% uncertainty in the calculated width of 
the ground-state molecular rotation and vibration distribution would contribute 6 x 1CT 3 eV 2 
to the budget for cr sys t (m^) [6]. 

Other sources of systematic uncertainty for KATRIN are more amenable to experimental 
control [6]. An electron gun behind the 10 11 Bq windowless, gaseous T 2 source will allow 
calibration of the experimental transmission function and of the energy loss experienced 
by electrons traveling through the source. The retarding potential of the KATRIN MAC- 
E filter will be independently monitored by the refurbished spectrometer from the Mainz 
experiment [35] and by a high-voltage divider with a demonstrated stability of 6.0 x 10 - ' 
per month [36]. Fluctuations in the column density of the source, which affect the scattering 
probability for electrons exiting the source, will be limited to the 0.1% level through control 
of the tritium injection rate, the pumping speed, and the vessel temperature; a temperature 
stability of 5 x 10 -5 per hour at 30 K has been demonstrated with a prototype system [37]. 
In addition to the primary component T 2 , it is expected that the KATRIN source will also 
contain DT and, to a lesser extent, HT. To achieve the desired stability of the column density 
and column activity, the isotopic purity of the source must be determinable to a relative 
precision of 2 x 10 -3 [6], and to this end the composition of the source gas will be monitored 
via Raman scattering in the tritium recirculation loop that feeds the source [38, 39]. 

Today, the beta decay of molecular tritium provides the most immediate path to im¬ 
proving the sensitivity to neutrino mass by direct, laboratory determination. Both the 
anticipated sensitivity of the KATRIN experiment now under construction and the develop¬ 
ment of new ideas motivate a careful evaluation of the 3 HeT + states excited in tritium beta 
decay. 


III. FORM OF THE BETA SPECTRUM 


The tritium decay process is accurately described by the Fermi theory of beta decay [40] . 
Tritium and helium-3 are mirror nuclei, so the nuclear matrix element M nuc is maximal. The 
transition is allowed, and the spectrum is not significantly modified by a shape factor depen¬ 
dent on the kinetic energy of the electron. Hence the shape of the beta decay spectrum is 
determined by the neutrino mass m„; electron mass m e ; total electron energy E e ; maximum 
energy of the electron, F max = Q — E^™ + m e \ and the energies 14 and probabilities P% asso¬ 
ciated with excitations of the daughter ion. The recoil energy E^™ consists of translational 
kinetic energy of the daughter ion. Since the discovery of neutrino oscillations shows there 
are three different neutrino eigenmasses (m„j), the full spectrum becomes an incoherent sum 
over individual spectra for mass index i — 1,2, 3, with intensities given by the squares of 
the neutrino mixing matrix elements (U e j) [41]. The resulting distribution of the electron 
energy E e is shown in Eq. 1, in which Gf is the Fermi weak-coupling constant, 9c is the 
Cabibbo angle, F(Z, E e ) is the Fermi function correcting for the interaction between the 
electron and the nucleus, and 0(A max — E e — 1% — m ui ) is a Heaviside step function ensuring 
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energy conservation [10]. Units are chosen where c — 1. 
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A number of small corrections to this basic spectral form have been identified over the years 
and have been summarized by Wilkinson [42] . At the time of his work, the effects he enu¬ 
merated were for the most part negligible, but as experimental precision has advanced, their 
significance has as well. Radiative corrections are the most important and have subsequently 
been re-examined [43]. A comprehensive and fully relativistic treatment of weak magnetism 
and induced terms may be found in Ref. [44]. 

Formally, Eq. 1 also contains inaccuracies in its treatment of rotational and vibrational 
molecular excitations. The mass of the nucleus is considered to be infinite in deriving the 
electron-neutrino phase space, and nuclear recoil is then treated separately in determining 
the molecular translation, rotation, and vibration in the final state. Electronic excitations 
represent energy unavailable to the outgoing leptons, and the modification to the phase 
space is appropriately captured by the appearance of 14 in expressions for the electron 
energy, ffowever, a correct treatment of rotational and vibrational excitations becomes 
ambiguous inasmuch as the appropriate recoil mass is not defined. In addition, the center- 
of-mass frame invoked for the decay described by Eq. 1 is not related in any simple way to 
the center of mass of an object more complex than an isolated atom. In a molecule, the 
atoms are always in motion, a source of Doppler broadening for the observed electron. These 
issues can be avoided by consideration in a relativistic formalism of the full three-body phase 
space populated in the decay. 

Because of the momentum imparted by the leptons to the recoil nucleus, the phase space 
is three-body rather than two-body everywhere except at the endpoint. While it is standard 
to neglect this effect, doing so introduces a small spectral distortion. More importantly, the 
three-body form permits a self-consistent treatment of recoil effects. The spectrum endpoint 
is given without ambiguity for any molecular system by conservation of the four-momentum 
for the full system. An exact three-body, fully relativistic calculation for the phase-space 
density has been given by Wu and Repko [45] (see also Masood et al. [46] and Simkovic et 
al. [44]): 
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with the following definitions: 
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We have here generalized Wn and Repko’s result by introducing multiple neutrino mass 
eigenstates • The mass M is the mass of the initial (final) atom or molecule, 

including associated atomic electrons and any excitation energy that may be present. The 
quantity A,, an experimentally useful fit parameter, is the ‘extrapolated endpoint energy’ 
that is obtained when the neutrino mass in the term in square brackets in Eq. 2 is set to zero. 
The quantity E eiraajK is the maximum energy of the electron for each neutrino eigenmass [47] . 
The electron-neutrino correlation modifies the spectrum at recoil order (~ m e /M ) [44] and 
is not included here. 

Both initial- and final-state excitations can now be introduced explicitly by indexing M 
and Mtf\ to become Mj and M^k, respectively. For each pair of initial and final states jk 
there is a corresponding Q-value, 


Qkj = Mj - M (f)k - m e (7) 

which is the kinetic energy released in the transition in the absence of neutrino mass. A 
special case is the atomic mass difference between the neutral atoms T (mass M 0 = A) and 
3 He (mass AI(j)o + m e = A') in their ground states, which we denote Qa- 

Qa = A- A'. (8) 

This corresponds to the Q-value for bound-state beta decay from ground state to ground 
state, the kinetic energy being delivered entirely to the neutrino and recoil. (The term 
“Q-value” without qualification is used inconsistently in the literature, sometimes meaning 
Q oo and sometimes Qa- For the atomic case, those quantities differ by the single-ionization 
energy of He, 24.59 eV.) 

In the general case the masses Mj and M(j\ k can be related to atomic masses by accounting 
for electron binding energies and for the possible presence of other atoms in the molecule: 

Mj = A s + A — bj (9) 

Af(/)fc = A s + A! — b(f) k — m e (10) 

Qkj — Qa bj -)- (ll) 

Here, the binding energies bj and b(f) k are the energies released in transforming an atomic 
mass to the species of the parent or daughter, and the atomic mass of the other, ‘spectator,’ 
nucleus in the molecule (if present) is denoted A s . For example, the binding of two neutral 
tritium atoms to form a neutral T 2 molecule in its ground state occurs with the release of 
b 0 = +4.59 eV. Figure 2 is a graphical summary of the relevant binding energies. 
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FIG. 2. Energy levels relevant to atomic and molecular tritium decay, patterned after Fig. 5 in 
Otten and Weinheimer [7]. The mass difference Qa is taken from Audi, Wapstra, and Thibault [48]. 
Dissociation energies are derived from calculations by Doss [49]; the ionization energy of 3 HeT + is 
from calculations by Kolos et al. [50]. The ionization energies for T [51] and for 3 He [52] are taken 
from recent compilations. 


The extrapolated endpoint energy Acan be expressed in terms of the corresponding 
Q-value: 

^ikj = Qkj + m e ~ 2 M■ + ^ me ) _ 2 M- ' (^) 

The extrapolated endpoint still has a dependence on neutrino mass, but it is completely 
negligible so the mass-eigenstate subscript i on A will be omitted henceforth. The recoil- 
order term is small, a few parts in 10 4 of Qk r Thus the extrapolated endpoint energy A k j 
for excited final states (< 100 eV) can be taken to be the ground-state quantity A 0 j minus 
the excitation energy. 

Weighting each transition by a matrix element Wkj for the transition connecting the 
specific initial state j to the final state k, the spectral density becomes 
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An expression for the matrix element W k j is given in Eq. 15 in Sec. IV. 

The maximum kinetic energy HA]" max ( fcj ) imparted to the recoil atom or molecule is the 
difference between the extrapolated endpoint energy and the available mass energy in the 
decay: 

■^rec, max ( kj) = + ^ m e) (14) 

A correct evaluation of the recoil energy is important because, as will be shown, the variance 
of the final-state distribution in the ground electronic state is directly proportional to it. 
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Tabic I summarizes the values of these parameters for several parent species, evalu¬ 
ated using the atomic mass difference Qa = 18591.3(10) eV given by Audi, Wapstra, and 
Thibault [48]. In ref. [53] a more recent measurement and a discussion of the experimental 
status of Qa are presented. 

TABLE I. Values in eV of the binding energies, Q-values, extrapolated endpoint energies, and 
maximum recoil translational energies for five tritium-containing parents. All of the quantities in 
the last three lines have the fractional uncertainty of Qa- 





Parent 



Vc^ u. cHILl 

x + 

T 

HT 

DT 

t 2 

bo 

-13.61 

0 

4.53 

4.57 

4.59 

b (f) 0 

-79.01 

-24.59 

-11.77 

-11.73 

-11.71 

Q 00 

18525.85 

18566.66 

18574.96 

18574.95 

18574.95 

A 0 o - rn e 

18522.44 

18563.25 

18572.40 

18572.91 

18573.24 

jrp kin 

■^reCjinaxfOO) 

3.402 

3.409 

2.557 

2.045 

1.705 


In particular, it may be seen from the table that the endpoint energy for HT falls about 
0.8 eV below that for T 2 , and the endpoint energy for DT is intermediate between the 
two. However, the same underlying kinematics produce a compensating energy shift in the 
final-state distribution, as described in Sec. V. 


IV. THEORY OF MOLECULAR TRITIUM BETA DECAY 

Molecular states are specified by electronic (n), vibrational (u), rotational (J), and az¬ 
imuthal (M) quantum numbers. For homonuclear molecules such as T 2 the total nuclear 
spin (/) is important in satisfying the Pauli exclusion principle. The T 2 nuclear spin can take 
on two values, 1 and 0; / = 1 corresponds to the triplet ortho state, and 1 = 0 corresponds 
to the singlet para state. The relevance of ortho and para states to the rotational quantum 
number and true molecular ground state is discussed in detail in Sec. IV D. 

The final states excited in molecular beta decay include translational, electronic, ro¬ 
tational and vibrational excitations. For the beta decay of an isolated tritium ion, only 
translational recoil is possible. For a neutral tritium atom, precisely calculable electronic 
excitations also occur. For a tritium molecule, rotational and vibrational excitations come 
into play and a theoretical treatment requires extensive computation. Even for a parent 
molecule as simple as T 2 , the electronic excited states of the daughter 3 HeT + molecule are 
complicated and unbound. Experimental advances, however, allow an important simplifica¬ 
tion: high statistics and excellent energy resolution will allow KATRIN to concentrate data 
taking within about 20 eV of the electron endpoint, a region in which electronic excitations 
play no role. Theoretical work can then focus on a precise calculation of the rotational and 
vibrational state distribution within the electronic ground state manifold. 

High-precision, ab initio calculations of the molecular excitations arising from T 2 beta 
decay have been performed [8, 9]. The calculations use the Born-Oppenheimer approxi¬ 
mation to factorize the molecular wavefunctions into electronic wavefunctions, vibrational 
wavefunctions and spherical harmonics dependent on the rotational and azimuthal quantum 
numbers. Hyperhne structure is neglected except where spin symmetry must be respected 
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in homonuclear systems. Corrections to the Born-Oppenheimer and other approximations 
have also been investigated and found to be small [54], 


A. Geminal-basis method 

Theoretical investigations of beta decay in T 2 date back to the pioneering work of Cantwell 
in 1956 [55]. Modern calculations are built on the theoretical framework of Kolos and Wol- 
nicwicz, who developed an adiabatic description of the hydrogen molecule in a basis of ex¬ 
plicitly correlated two-electron wavefunctions in 1964 [56]. This basis is sometimes described 
as geminal because it treats the electrons as a pair rather than as independent particles. 
Development of the geminal basis for the hydrogen molecule led to early calculations of the 
molecular effects in the decay of HT [57]. In a further refinement of the basis, Kolos et 
al. [50] investigated optimal parameter values. The most recent calculations rely on those 
results with minor additional refinements [8]. 

As Jonsell, Saenz, and Froelich [54] show, the transition matrix element related to the 
final-state 3 HeT + excitation k = (w(/), </(/), from an initial T 2 state j = (v, J, M, n) 

may be written, 


\W ki (K)f 


Hf) 




R) 


S n (R)e^^ JM (R)d^R 


(15) 


In this expression, x an d £ are the rotational-vibrational wave functions of the 3 HeT + and 
T 2 molecules, respectively, and S n (R ) is an electronic overlap integral. The exponential of 
the dot product of the recoil momentum K and the nuclear separation R is a consequence 
of the recoil motion of the daughter He nucleus. 

The reduced mass of the daughter molecule enters into the radial Schrodinger equation, 
which must be solved in order to compute the rotational and vibrational energy levels. 
There is some ambiguity in the definition of this quantity, which depends on whether and 
how the masses of the two bound electrons are included. Coxon and Hajigeorgiu [58], 
comparing predicted energy levels to spectroscopic measurements (Sec. VIA 1), achieved the 
best agreement with an effective reduced mass that assumes one electron belongs strictly to 
the He nucleus, with the second electron distributed evenly between the H and He nuclei. 
Doss et al. [9], confirming this result, introduced the effective reduced mass to the calculation 
of the final-state distribution, but noted that the change was insignificant at the 0.1-eV level 
of foreseeable T 2 -based neutrino-mass measurements. 

Fig. 3 shows the spectrum of final-state molecular excitations from the beta decay of 
T 2 (J = 0 initial state) as published by Saenz et al. in 2000 [8], compared with the 1985 
calculation by Fackler et al. [59]. The electronic ground state appears as a large peak 
centered at 1.7 eV excitation energy (0.2 eV binding energy), broadened by the rotational 
and vibrational excitations. The higher electronic states also suffer broadening as shown. 
For a detailed description of the differences between the Saenz and Fackler calculations see 
reference [8] . The more recent results of Doss et al. [9] were not published in tabular form 
but a subset of the tables was provided to the KATRIN collaboration courtesy of Doss. 
Reference [49] compares the Doss et al. [9] results and the Saenz et al. [8] results. The 
differences are negligible for the ground state but noticeable in the electronic continuum, 
particularly above 45 eV of excitation energy (see Sec. IV C). 
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FIG. 3. Molecular spectrum excited in the beta decay of T 2 (J = 0) as calculated by Saenz et al. [8] 
(solid curve, red online) and by Fackler et al. [59] (dotted curve, blue online). For the purposes of 
display and comparison, discrete states in the latter spectrum have been given a Gaussian profile 
with a standard deviation of 3 eV. 


Unfortunately, in the geminal basis the convergence of the calculations depends on the 
explicit choice of basis functions and in-depth study has revealed that adding even one 
basis function can dramatically change the contributions of other functions [60]. Significant 
optimization work was done to select the correct basis functions for T 2 and 3 HeT + and 
provide reliable results [50] . I 11 lieu of explicitly computing uncertainties, which is impractical 
due to the volatility of the basis, published calculations typically include the deviation from 
1 of the cumulative probability function at the highest excitation energy. However, this 
single number, while informative, is an insufficient gauge of accuracy. Despite the linear 
dependences and instability of the geminal basis, it has been used to produce essentially all 
final-state distribution calculations of the ground-state manifold since its publication [8, 9]. 


B. Configuration-interaction method 

The configuration-interaction (Cl) method presents an alternative approach to model¬ 
ing two-electron, diatomic molecules such as T 2 and 3 HeT + within the Born-Oppenheimer 
approximation. In the Cl technique, two-electron configurations are constructed as the 
products of pairs of solutions to the single-electron Schrodinger equation (denoted orbitals). 
Superpositions of these configurations are then used to build wavefunctions and make calcu¬ 
lations. As the simplest two-electron heteronuclear molecule, HeH + was an early test bed for 
the method (see, e.g., [61-63]). I 11 the 1980s, parallel to the refinement of the geminal-basis 
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method, the Cl method was applied to the FSD following beta decay in T 2 . Fackler et al. [64] 
performed a preliminary study of decays to the first five electronic states of 3 HeT + ; Martin 
and Cohen [65] used a more flexible basis set of Cartesian Gaussian orbitals to study the first 
50 eV of the electronic continuum. (See Sec. IV C.) Without the benefit of modern compu¬ 
tation, however, such early treatments were neither complete nor precise enough to address 
the final-state spectrum in the region of interest for modern tritium-based neutrino-mass 
experiments such as KATRIN. 

More recently, Vanne and Saenz [66] have developed a Cl approach, based on an un¬ 
derlying B-spline basis set and carried out in an elliptical box, that shows promise for 
neutrino-mass experiments. This method avoids the linear dependences that tend to arise 
in numerical calculations with the geminal-basis method, allowing application to larger inter- 
nuclear distances R as well as the use of larger basis sets. Adding individual basis functions 
does not introduce artificial resonances. The discretization provided by the elliptical box al¬ 
lows the electronic continuum to be discretized as well, permitting the consideration of both 
bound and continuum states within the same basis set. Since all configurations are expressed 
in terms of one-electron wavefunctions, however, two-electron correlations are treated less 
accurately than in the geminal-basis method, especially if the configuration set is small. 

Vanne and Saenz have compared their B-splinc-based Cl treatment of HeH + photoion¬ 
ization [66] against one using the standard geminal basis [67]. The first resonance in the 
X 1 ^ —y photoionization cross section, at about 16 eV, is shifted about 0.5 eV higher 
in the Cl results, likely due to the difference in treating two-electron correlations. The two 
approaches predict the same amplitude for this resonance and give good agreement for other 
features of the spectrum. 

The application of this method to tritium beta decay is a work in progress [68]. Once 
sufficiently complete configuration sets are calculated for T 2 and for 3 HeT + , the electronic 
overlap integrals S n (R) can be computed. Transition probabilities may then be determined 
using Eq. 15. 


C. Electronic continuum 

The energy window for the KATRIN neutrino-mass measurement is narrow enough that 
related FSD calculations can focus on the 3 HeT + electronic ground state. However, it has 
been suggested that a measurement of the tritium beta spectrum over a wider energy range 
could be used to search for sterile neutrinos with mass on the eV scale [69] or even on the 
keV scale [70]. If the acquisition window extends more than about 40 eV below the beta 
endpoint, the analysis must account for the electronic continuum portion of the FSD. Table II 
gives a brief overview of the variety of methods that have been applied to the problem. In 
addition to their differences in general approach, the available calculations differ in baseline 
assumptions. Early calculations often used the clamped-nuclei approximation rather than 
explicitly accounting for nuclear motion that broadens resonances. Assumptions about the 
localization of resonances can introduce errors at higher excitation energies [71]. Variation 
of the internuclear distance shifts the overall probability distribution but can also change the 
relative intensities of the electronic resonances [72], A significant simplification is possible 
at excitation energies above ~ 200 eV, a region in which the fast-moving ejected electron 
sees the 3 He ++ ion as equivalent to a bare He nucleus. The high-excitation-energy tail of the 
FSD can then be described with a spectrum adapted from the decay of atomic tritium [72], 

The calculated percentage of tritium decays that populate the electronic continuum is 
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TABLE II. Selected calculations of the probability P C ont of populating the electronic continuum of 
3 HeT+ in T 2 beta decay. The integration range differs between the calculations, and the bounds 
are specified as excitation energies above the 3 HeT“ t ~ ground state. 


Method 

Reference 

Pcont 

Integration Range 

Complex scaling 

Froelich et al. (1993) [73] 

12 . 77 % 

45 - 90 eV 

Stieltjes imaging 

Martin and Cohen (1985) [65] 

13.42% 

45 - 94 eV 

Stabilization 

Fackler et al. (1985) [59] 

14.2% 

45 - 200 eV 

R-matrix 

Doss and Tennyson (2008) [71] 

13.66% 

ca. 40 - 240 eV* 


*Lower integration bound is not explicitly given. 


relatively consistent despite dramatic differences in the integration range, reflecting the fact 
that this region of the spectrum is dominated by a few autoionizing states near the ionization 
threshold. However, comparisons between different calculations, performed e.g. in Ref. [71] 
and [73], show significant discrepancies in the detailed structure of this part of the spectrum. 
For a sterile-neutrino search, knowledge of the integrated probability P con t is not sufficient. If 
not properly accounted for, small structures in the FSD at high excitation energies could lead 
to errors in interpretation, especially when small mixing angles are considered. Sensitivity 
calculations for such a search must be guided by theoretical studies of this region of the FSD 
spectrum. 


D. Molecular forms of tritium 

The tritium-containing hydrogen isotopologs (HT, DT and T 2 ) have different reduced 
masses and thus different excitation spectra. While the overall structure of the final-state 
spectrum remains qualitatively the same across isotopologs, the vibrational energy levels 
are shifted and the probability of a transition to any specific rotational-vibrational state 
changes. For example, the electronic excitations in 3 HeH + are shifted ~1 eV lower than the 
corresponding excitations in 3 HeT + [54], As shown in Table I, however, the difference in 
recoil mass also changes the extrapolated endpoint, canceling the change in the beta energy 
to first order [6]. 

In addition to differences in reduced mass, nuclear spin and symmetry considerations play 
an important role in determining the allowed angular-momentum states of the homonuclear 
T 2 molecule but do not apply to the heteronuclear DT and HT molecules. In accordance 
with Fermi statistics, the overall T 2 wavefunction must be antisymmetric under exchange of 
the tritium nuclei. The electronic, rotational, and vibrational wavefunctions of the molecule 
are inherently symmetric. Thus the spin-symmetric ortho state must be matched with an 
antisymmetric spatial wavefunction corresponding to odd J. The spin-antisymmetric para 
state must be matched with a symmetric spatial wavefunction corresponding to even J. 
Hence the ground state of the molecule is the para state with J = 0. 

In thermal equilibrium the partition function of rotational states ( J) in T 2 may be written, 

OO 

^equii = - (-1) J ](2J+ l)e- J ( J+1 ) R2 / 2 -' fcflT , (16) 

j =0 

to first order. Here the first factor is the spin statistical weight for ortho (odd J ) or para 
(even J) in the case of a homonuclear molecule, when total antisymmetry must be enforced, 
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and ksT is the thermal energy. The moment of inertia, , is related to the energy of the 
first excited state, E J=i , 




h 2 

2E J= i 


(17) 


Since Ej = i = 0.00497 eV [49] is small compared to hsT at room temperature, the ortho¬ 
para ratio of a thermally equilibrated source at room temperature is essentially the ratio 
of the spin statistical weights, 3:1 [54], Rather than the ortho-para ratio, the state of a 
molecular hydrogen source is typically characterized in terms of the parameter A quantifying 
the fraction of the source that is in the ortho state. 

The ortho-para transition requires a simultaneous change in the spin and rotational quan¬ 
tum numbers, making the ortho state metastable. Thus transitions to lower rotational states 
are dominated by intrinsically slow quadrupole transitions. For this reason, unless specific 
steps are taken to ensure it, thermal equilibrium of the rotational states of T 2 cannot be 
guaranteed. Thermalization of the spin degrees of freedom in a homonuclear hydrogen source 
is a slow, exothermic process, and uncertainty arises from the use of sources that are not in 
thermal equilibrium and that contain a mixture of states. 

Previous studies of molecular hydrogen have focused on the ortho-para ratio alone as 
the determining factor in the rotational-state distribution, a reasonable assumption for light 
isotopologs. However, for T 2 above cryogenic temperatures, states higher than J = 1 have 
significant populations and the evolution of the full rotational-state distribution must be 
considered. Spontaneous quadrupole transitions are extremely slow, on the order of ICR 7 s -1 
in free space [74], and transitions will be dominated by collisions with other tritium molecules 
and the walls. The rate of these processes depends on the detailed design of the gas system 
and must be carefully modeled to determine the rotational-state distribution of the source. 


V. CONCEPTUAL MODEL OF THE ROTATIONAL-VIBRATIONAL SPEC¬ 
TRUM 

As we have seen in Sec. IV, a precise treatment of the molecular final-state spectrum 
requires an extensive theoretical framework. However, as experimental sensitivity has ad¬ 
vanced, dependence on the highly excited states has diminished. The width of the ground- 
state manifold now sets the fundamental limit on the sensitivity of experiments using T 2 . 
With the intention of gaining some insight into the physical origin of the width of this 
manifold we have developed a simplified treatment, based on kinematic considerations and 
the approximation of the molecule as a simple harmonic oscillator. It reproduces several 
features of the precisely calculated spectrum while clarifying the underlying physics. 

Qualitatively, the beta spectrum is influenced in two distinct ways by the molecular 
structure. The rotational, vibrational and translational motions of the parent T 2 molecule 
lead to modulation of the energy of the detected beta electron. Some motions are essen¬ 
tially thermal in origin and contribute a Doppler shift in the laboratory electron energy. 
Classically, each degree of freedom contains on average \ksT of energy, and the atomic ve¬ 
locity adds vectorially to the electron velocity. Nevertheless, as we shall see, it is a uniquely 
quantum-mechanical effect, zero-point motion, that in fact dominates the spectrum at low 
temperatures. 

In the following, our interest is in the rotational and vibrational degrees of freedom in 
the electronic ground state. We begin by examining the purely kinematic constraints on 
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the recoil momentum p = KK. We then, in a semi classical approach, combine the initial 
momentum of the decaying T nucleus in the parent molecule with the momentum delivered 
by lepton recoil in order to find the momentum spectrum of the daughter 3 He. Applying 
kinematic constraints, the momentum spectrum is expressed in terms of the corresponding 
translational and excitation energies of the recoil molecular ion 3 HeT + or 3 HeH + , for the 
parents T 2 and HT, respectively. 


A. Recoil momentum 

The three-momentum imparted to the molecular system by the beta decay has a magni¬ 
tude 


V = |Pe + p„| 

P 2 = E 2 e - m 2 e + (£ max ~E e - E^) 2 - ml + 2 E e (E max - E e - A r k j c n )/3^ cos 9 ev (18) 

where 9 ev is the angle between the electron and the neutrino momenta, and /3 and j3 v are, 
respectively, the electron speed and neutrino speed relative to the speed of light. It is 
sufficient for the present purpose to neglect neutrino mass and also the kinetic energy of the 
recoil as it contributes corrections of order m e /M ~ 10 -4 to the square of the recoil 
momentum. 

The electron-neutrino correlation term may be written [75] 


1 + a ev 


Pe ■ Pu 

E e E u 


1 + a ev /3 cos 9 ev . 


(19) 


Using for a eu the value measured for the free neutron, a ev = 0.105(6) [76], and noting that 
the electron velocity (3 < 0.26, one sees that the electron-neutrino correlation is very weak 
in tritium decay. The recoil-energy envelope for the decay of an isolated tritium nucleus is 
shown in Fig. 4. 

Although the recoil momentum is given immediately from the lepton momentum via 
momentum conservation, determining the recoil energy requires knowledge of the recoil 
mass. In the case of an isolated T atom, shown in Fig. 4, the calculation is unambiguous, 
but for a T 2 molecule it is not. For a very tightly bound system with no accessible internal 
degrees of freedom the mass would be the total mass (6 u), and for a very weakly bound 
one it would be 3 u. Without further information, the recoil energy can be bounded above 
and below by kinematics and at these limits is entirely translational kinetic energy. At the 
endpoint of the beta spectrum, 


1.705 < Ely < 3.410 eV. 


( 20 ) 


The 3 HeT + ion has a spectrum of rotational and vibrational excitations that are one or 
two orders of magnitude smaller than the recoil energy, less like the strongly bound picture 
and more like the weakly bound one. Some insight into the behavior of this system can be 
gained by considering first a purely classical T 2 molecule at 0 K, such that both atoms are 
bound together but at rest. If the molecule remains bound after beta decay, conservation 
of linear momentum requires that 1.705 eV must be in the form of translational kinetic 
energy, leaving only 1.705 eV available for internal excitations. The binding energy of the 
final-state molecular ion 3 HeT + is 1.897 eV [49], and, since this is greater than the available 
excitation energy, the 3 HeT + must remain bound in this classical picture with no thermal 
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FIG. 4. Recoil kinetic energy imparted to a 3 He daughter by the beta decay of an isolated tritium 
nucleus at rest. The upper boundary of the envelope corresponds to 9 ev = 0 and the lower one to 
9ei/ = If ■ 




FIG. 5. Distributions of excitation energy in the ground-state rotational and vibrational manifold 
of 3 HeT + (left) and 3 HeH + (right) as calculated by Saenz et al. [8]. The expected value for the 
excitation energy in each case, based on kinematic considerations, is indicated by a vertical line. 
An excitation energy of 0 corresponds to a binding energy of 1.897 eV [49]. 


motion. Then the final state consists of a mass-6 ion with a translational kinetic energy of 
1.705 eV and rotational and vibrational excitations totaling 1.705 eV. How the excitation 
energy is apportioned between rotational and vibrational excitations depends (classically) on 
the relative orientation of the axis connecting the atoms to the lepton momentum direction, 
but the total excitation energy is always 1.705 eV. 

The equivalent relationship for the HT parent molecule is 

2.557 < < 3.410 eV, (21) 
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the total internal excitation of the 3 HeH + is 0.85 eV, and the translational kinetic energy 
is 2.557 eV. We compare these expectations with the calculations of Saenz, Jonsell, and 
Froelich [8] in Fig. 5. 

The centroids of the theoretical distributions agree with our expectation but the distri¬ 
butions are not delta functions. Broadening is caused by the fact that atoms in the parent 
molecule are always in motion due to thermal and quantum effects, which smears the final- 
state momentum of the 3 He and the momentum of the outgoing leptons. The calculations 
of Saenz et al. were carried out in the center of mass for T 2 and HT gas at 30 K; we shall 
show that, at low temperatures, the chief mechanism for broadening is zero-point motion in 
the parent molecule. 


B. Spectrum of the electronic ground state 

A diatomic molecule at low excitation may be described as a one-dimensional harmonic 
oscillator: 


E v = (v + 1/2 )hhJ c + a(vhu > c ) 2 ; v = 0,1, 2,... (22) 

(23) 

where k is the force constant for displacements from the equilibrium internuclear separation, 
and /i is the reduced mass. A small anharmonic term with coefficient a is included. By fitting 
the four lowest vibrational states of the H 2 molecule [77] one folds hco c = 0.5320(5) eV and 
a = —0.0537(8) eVW The corresponding value of hcu c for T 2 is then 0.3075 eV, much larger 
than ksT at 30 K (0.003 eV), and also larger than typical rotational excitations (0.005 eV). 
In the vibrational ground state, the zero-point motion has an equivalent temperature of 
about 0.15 eV, or ~ 1600 K, and dominates the line broadening. The zero-point energy is 

(24) 

When beta decay occurs, the lepton recoil momentum p adds vectorially to the instan¬ 
taneous momentum p T of the decaying tritium nucleus of mass m within its molecule: 




p f = p + p T . 

The mean kinetic energy of the decaying tritium nucleus is 


M> = l -R 

2/i 


/i = 


2 ~ zp ’ 
m q m 


m s + m 

and the standard deviation of the excitation energy E exc of the recoil ion is then 


= -JUpP 

mV 3 


pZ 


' 2 m l 3m zp 


2 fi 


(25) 

(26) 

(27) 

(28) 
(29) 
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FIG. 6. Calculated recoil excitation energy spectra from zero-point motion in the parent molecule 
(dotted curves, blue online), compared to the final-state distributions calculated by Saenz et al. [8] 
at 30 K (solid curves, red online). The curves from zero-point motion are parameter-free except 
for normalization, and have the standard deviations indicated in Table III. An excitation energy 
of 0 corresponds to a binding energy of 1.897 eV. 


where m is the mass of the decaying tritium nucleus, and m s is the nuclear mass of the 
‘spectator’ nucleus in the molecule. For the present purposes we ignore the difference between 
the nuclear masses of T and 3 He. 

Inserting for E zp the relevant zero-point energies for T 2 and HT, the predicted distri¬ 
butions of recoil excitation energy are compared with the calculated spectra of Saenz et 
al. [8] in Fig. 6. The good agreement (4%; see Table III) underscores the fact that the gross 
features of the final-state distribution really arise from the initial state, i.e. it is mainly the 
zero-point motion of the tritium atom in its molecule that broadens what would otherwise 
be a line feature. The broadening occurs even for the ground-state molecule at absolute 
zero and is irreducible. Final-state effects assert their presence only through the density of 
available states in the 3 HeT + and 3 HeH + ions, which modulates the continuous distribution. 
That modulation may be calculated by evaluating the overlap integral between the final- 
state wave functions and the momentum projection operator acting on the initial state as 
given above. 


TABLE III. Root-mean-square widths in eV of the ground-state manifold from the exact calculation 
of Saenz et al. [8] for initial state J = 0, and derived from the semiclassical treatment based on 
the zero-point motion of the parent molecule. 


Method 

t 2 

HT 

Saenz et al. [8] 

0.436 

0.379 

Zero-point motion 

0.420 

0.354 


Including the smearing effect of zero-point motion, the line feature is broadened to a 
sufficient extent that a large fraction of the distribution lies above the dissociation threshold, 
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1.897 eV in 3 HeT + . Jonscll et al. [54] find that while the intensity of the ground state 
transitions near the T 2 beta endpoint amount to 57% of decays, 18% absolute lies above 
the dissociation threshold. For HT only 1.5% absolute is above the dissociation threshold 
(see Fig. 6). Not all such excited states will necessarily dissociate, however, because of the 
angular momentum barrier for states with high J. Those states may be sufficiently long-lived 
to decay radiatively. 

The T 2 vibrational energy interval of 0.308 eV is almost two orders of magnitude larger 
than the excitation energy Ej = \ = 0.00497 eV of the lowest-lying ortho state (Sec. IV D); 
the zero-point motion is thus the dominant contribution to the final-state width. If the 
parent molecule is in an initial state with angular momentum J, the root-mean-square 
width becomes 


^"(-^exc) 


, 2 a 2 m 2 e J(J + 1) \ 
2 m \3 m zp 3Rim / 


(30) 


where a is the fine structure constant and 7% is th e equilibrium internuclear separation in 
a.u. (1 a.u. = h/m e a). The variances of the excited-state distributions for T 2 , DT, and 
HT for states up to J = 10 are given in Table IV, and a graphical comparison with the 
calculations of Doss [49] up to J = 3 is displayed in Sec. VII. 


TABLE IV. Root-mean-square widths in eV of the ground-state manifold of the daughter molecule 
from the semiclassical treatment based on the zero-point motion of the parent molecule with the 
inclusion of rotation. 


(v,J) 

t 2 

DT 

HT 

(0,0) 

0.4197 

0.3972 

0.3537 

(0,1) 

0.4331 

0.4113 

0.3694 

(0,2) 

0.4586 

0.4381 

0.3991 

(0,3) 

0.4944 

0.4755 

0.4398 

(0,4) 

0.5385 

0.5212 

0.4888 

(0,5) 

0.5890 

0.5732 

0.5439 

(0,6) 

0.6443 

0.6299 

0.6034 

(0,7) 

0.7035 

0.6903 

0.6662 

(0,8) 

0.7654 

0.7533 

0.7313 

(0,9) 

0.8297 

0.8185 

0.7983 

(0,10) 

0.8956 

0.8853 

0.8667 


The objective in a tritium beta-decay experiment is measurement of the neutrino mass 
via a detailed study of the shape of the electron spectrum near the endpoint. Energy 
conservation assures a connection between the molecular final state and the electron energy. 
The modification can be directly derived and has a particularly appealing and simple form. 

If the tritium atom has a velocity fir in the center of mass at the instant the decay takes 
place, the foregoing considerations of zero-point motion in the molecule give 

(ft) = ^ ms 

' ' 3m m s + m 

a(E e ) = E e /3yf(ft). (31) 

This result is identical to Eq. 29, the previously derived width for the excitation of the recoil. 
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C. Recoil energy spectra in dissociation 


The theory of molecular beta decay can also be used to predict the energy of the ions 
produced in the decay. A measurement of the ion energy spectra would be helpful in assessing 
our understanding of the underlying decay. As Sec. VIB 4 discusses in detail, theory predicts 
that approximately half of the decays of T 2 and HT lead to dissociative states [54], whereas 
experimental data indicate that more than 90% of the transitions lead to bound molecular 
ions [78, 79]. While there are several plausible experimental and theoretical explanations for 
this discrepancy, the disagreement motivates an examination of the dissociation-fragment 
spectrum that would be predicted by theory. A future experiment may be able to make a 
measurement of this spectrum, providing a new test of the theory. 

We examine the six dominant electronic configurations in the Born-Oppenhcimer approx¬ 
imation as given by Jonsell et al. [54], These configurations account for 84% of the intensity, 
with the remaining 16% coming from the electronic continuum. In the ground-state man¬ 
ifold there is a potential minimum that leads to binding of the 3 HeT + by almost 2 eV; all 
electronic excited states are monotonically repulsive with the exception of the first excited 
state, which has a shallow minimum far outside the Franck-Condon region. Rotational and 
vibrational states in the electronic-ground-state manifold are quasibound because of the 
potential minimum augmented by an angular momentum barrier. For this analysis we con¬ 
sider these quasibound states to be rotational and vibrational states of a bound (mass-6) 
ion which dissociates by tunneling through the barrier, analogous to fission. Conversely, ow¬ 
ing to the absence of a binding potential, molecular motion in the electronic excited states 
corresponds more closely to the unbound scenario in which all the lepton momentum is de¬ 
livered to a mass-3 recoil ion. In this case the two fragments gain additional kinetic energy 
at dissociation by converting the repulsive potential energy of the excited molecular state 
at the Franck-Condon spatial separation. The necessary data for the latter calculation can 
be found in Fig. 1 of Ref. [54], 

The laboratory energies of the dissociation fragments from the quasibound ion can be 
calculated from kinematics. The laboratory kinetic energy for a fragment of mass m; 

is uniformly distributed in the interval 



(32) 


for A exc > Eb and rriiE^™ > rrij(E exc — Eb), where rrij is the mass of the other fragment and 
Eb is the binding energy of the molecular ion. It may be seen from this that the dissociation 
fragments from the quasibound states do not have translational energies significantly greater 
than that of the mass-6 ion, E^. 

Decays populating the electronic excited states produce recoil fragments, at least one 
of which is itself in an electronic excited state. Applying the Franck-Condon principle, the 
electronic excitation energy of the system before dissociation is evaluated at the internuclear 
separation of the T 2 , HT, or DT molecule in its ground state, 1.40 a.u. for all three species 
[9]. Tables V and VI list the relevant properties for the six dominant electronic states and 
the kinetic energies of the recoil fragments. 

The total kinetic energy available to the dissociation fragments by conversion of the 
interatomic potential in the five excited states is confined to a rather small range between 
10 and 17 eV. An additional amount of kinetic energy E [ trans ) = p 2 /2m is contributed to the 
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TABLE V. Structure of excited states and kinetic energies of dissociation fragments for the decay 
of T 2 . The probabilities, which are valid in the sudden approximation, are taken from [9] for the 
case Ji = 0, and are very similar for R = 1,2,3. The total probability calculated for these six 
states is 84.2%. 


State Asymptotic 

structure 

Excitation 
Energy (eV) 
R = 00 

Excitation 
Energy (eV) 
R = 1.40 a.u. 

Total Kinetic 
Energy (eV) 

E( 3 He) 

eV 

E(T) Branch 
eV (%) 

1 

3 He(ls 2 ) + T+ 

< 0 

< 1.897 

0 

0 

0 

39.0 


3 He(ls 2 ) + T+ 

> 0 

> 1.897 

Eq. 32 



18.4 

2 

3 He+(ls) + T(ls) 

10.981 

24.587 

13.606 

6.8 + 3.477 

6.8 

17.4 

3 

3 He(ls2s) + T+ 

20.5 

31.390 

10.890 

5.4 + 3.47/ 

5.4 

7.8 

4 

3 He+(ls) + T(2s+2p) 

21.186 

36.152 

14.966 

7.5 + 3.47/ 

7.5 

0.8 

5 

3 He+(ls) + T(2s-2p) 

21.186 

36.833 

15.647 

7.8 + 3.4/7 

7.8 

0.01 

6 

3 He(ls2p) + T+ 

21.0 

37.513 

16.513 

8.3 + 3.47/ 

8.3 

0.9 


TABLE VI. As Table V, for the decay of HT. The total probability calculated for these six states 


is 83.8%. 


State Asymptotic 

structure 

Excitation Excitation 
Energy (eV) Energy (eV) 
R = 00 R = 1.40 a.u. 

Total Kinetic 
Energy (eV) 

E( 3 He) 

eV 

E(H) Branch 
eV (%) 

1 

3 He(ls 2 ) + H+ 

< 0 

< 1.897 

0 

0 

0 

55.4 


3 He(ls 2 ) + H+ 

> 0 

> 1.897 

Eq. 32 



1.5 

2 

3 He+(ls) + H(ls) 

10.981 

24.587 

13.606 

3.4 + 3.47/ 

10.2 

17.4 

3 

3 He(ls2s) + H+ 

20.5 

31.390 

10.890 

2.7 + 3.47/ 

8.2 

7.8 

4 

3 He+(ls) + H(2s+2p) 

21.186 

36.152 

14.966 

3.7 + 3.47/ 

11.2 

0.8 

5 

3 He + (ls) + H(2s-2p) 

21.186 

36.833 

15.647 

3.9 + 3.47/ 

11.7 

0.01 

6 

3 He(ls2p) + H+ 

21.0 

37.513 

16.513 

4.1 + 3.4t/ 

12.4 

0.9 


recoil of the beta-decay daughter by the lepton momentum. We therefore define and use in 
the Tables a parameter 0 < // < 1 that is the fraction of the maximum lepton momentum 
squared. The other nucleus is a spectator and receives only the kinetic energy obtained from 
conversion of potential energy. The maximum energy imparted to a mass-3 daughter recoil 
is then about 12 eV for T 2 and 7 for HT. The He lines will be broadened by the zero-point 
motion as described in Sec. VB, and all lines will be broadened by the steep gradient of the 
interatomic potential in the Franck-Condon region. Moreover, in an experiment the total 
lepton recoil momentum is not directly measurable; only the electron momentum is. This 
introduces a range of values of 77 at each energy, as may be seen in Fig. 4 and Eq. 18. A 
detailed calculation of the line widths is beyond the scope of this work. 

The combination of the branching ratio to the bound molecular ion and the ion energy 
spectra provides a complete picture of the decay process. Measuring the branching ratio 
and kinematics has the potential to improve our understanding of the efficacy of our current 
model of molecular beta decay. 
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VI. TESTS OF TRITIUM FINAL-STATE CALCULATIONS 


The sub-eV energy scales of the rotational and vibrational excitations and the unknown 
time scales for further evolution of the final-state products make direct measurement of the 
final-state distribution, and particularly those aspects that are reflected in the corresponding 
lepton momentum, all but infeasible. Of particular concern are detector energy resolution 
and translational Doppler broadening of the distribution in a real experiment. The difficulty 
of a direct measurement has led to a variety of stratagems for indirect verification of the 
theory. In this section we discuss available data from spectroscopy, photodissociation, and 
mass spectrometry. 


A. Studies of the HeH + molecule 

1. Rotational and vibrational level transitions 

Determining the distribution of 3 HeT + final states populated by beta decay requires cal¬ 
culating the energy levels of T 2 and of 3 HeT + . If the same theoretical framework is also 
applied to calculating the spectra of molecules with other isotopes of He and H, predicted 
transition energies can be compared against a large number of transition lines measured 
with high-precision spectroscopic techniques ranging from glow discharge to absorption spec¬ 
troscopy to Raman spectroscopy. Such a comparison, of course, cannot test the probability 
of populating each 3 HeT + state after beta decay, but as we saw in Sec. IV A it has provided 
valuable input to modern theoretical calculations. 

Doss [49] calculated transition energies between rotational and vibrational levels in the 
electronic ground state for three tritium-containing parent molecules and for two daugh¬ 
ter molecular ions and compared them to published spectroscopic data. For 21 measured 
transitions in T 2 [80], seven in DT [80], and 12 in HT [80, 81], ranging between 120 and 
3775 cm -1 , the theoretical values always agreed within 1 cm -1 with a maximum fractional 
deviation of 0.1%. For 16 transitions in 3 HeH + and 10 in 3 HeD + [82-84], ranging from 71 to 
3317 cm -1 , the agreement is still better, within 0.05%. However, there do exist experimen¬ 
tally measured transition energies for which no geminal-basis predictions are reported: two 
rotational-vibrational Q 1 transitions in T 2 [85] and two in DT [86], three purely rotational 
transitions in the vibrational ground state of HT [86] , and 12 transitions in hot vibrational 
bands of HT [81] that fall well outside the energy range of the other measured transitions. 

In an earlier calculation in the standard geminal basis, Jonsell et al. [54] predicted transi¬ 
tion energies ranging from 598 to 3157 cm -1 in helium hydride molecular ions containing the 
more common isotope 4 He, allowing validation against a much broader catalog of spectro¬ 
scopically measured transitions. Five observed transition energies in 4 HeD + [82] and sixty- 
two in 4 HeH + [87-91] agree with these predictions to within 0.04%. The measured widths 
of seventeen predissociative resonances in 4 HeH + , 3 HeH + , 4 HeD + , and 3 HeD + [82, 91] differ 
from the predicted values by up to an order of magnitude, but the specific machinery for 
calculating these widths is not used to determine the final-state distribution for neutrino- 
mass measurements [54], No predictions are reported in the geminal basis for 46 additional 
observed transitions in low-lying vibrational bands of 4 HeD + [83, 84, 92, 93], or for 36 similar 
transitions in 4 HeH + [83, 84, 92, 94, 95]. 

Despite this great investment of experimental effort, only partial, fragmentary spectra 
have been measured for these seven molecules. Nonetheless, Coxon and Hajigeorgiu [58] 
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were able to use these data to construct a fitted Born-Oppenheimer potential for the generic 
molecular helium hydride ion HeH + , and compare it to an ab initio potential obtained from 
an older geminal basis with adiabatic corrections from Bishop and Cheung [96]. The two 
potentials differ by up to 2 cm -1 when the nuclei are close together but are in excellent 
agreement for internuclear distances R > 8 a.u.; the dissociation energies differ by only 
0.27 cm -1 [58]. No such comparison has yet been performed for the ab initio potential based 
on the most recent geminal basis. 

While theoretical predictions for all measured transition energies would be useful, the 
excellent agreement obtained over 133 transition energies in seven diatomic molecules sug¬ 
gests that the rotational and vibrational energy levels of the electronic ground states are 
well reproduced in the geminal basis. 


2. Photodissociation of^HeH + 

The photodissociation spectrum of 4 HeH + may be derived from a sufficiently complete 
theoretical description of the molecule. Since all electronic excited states of this molecule are 
dissociative in the Franck-Condon region, one can construct the photodissociation cross sec¬ 
tion as a function of energy by calculating dipole transitions between the electronic ground 
state and the electronic excited states. The result depends on the orientation of the inter¬ 
nuclear axis relative to the photon polarization vector; the parallel and perpendicular cases 
must be treated separately. Several other theoretical models ( e.g. [97, 98]) have been em¬ 
ployed to study the photodissociation problem, but have not been applied to neutrino-mass 
measurements. 

The process has been probed experimentally with 38.7-eV (32-nm) photons at the Free- 
electron LASer in Hamburg (FLASH). The initial measurement [99] determined the cross 
section to the He + H + channel, and was not able to define the initial distribution of vi¬ 
brational states in 4 HeH + . The second FLASH measurement [100] incorporated several 
experimental upgrades to provide additional tests. The 4 HeH + beam could optionally be 
routed through a linear electrostatic ion trap and cooled to the v = 0 vibrational ground state 
before being extracted to the interaction region. An improved detection setup, combined 
with a positive potential across the ion-photon interaction region, allowed the measurement 
of the branching ratio to the 4 He + H + and 4 He + + H channels. In both experiments, the 
distribution of the initial internuclear axis orientations was assumed to be isotropic. 

Beginning with the same geminal basis set as that used for standard neutrino-mass- 
relevant calculations, Saenz computed the total photoabsorption cross section assuming that 
the molecule begins with v = 0 and is oriented parallel to the photon held [67]. Dumitriu 
and Saenz later performed a more detailed calculation in the Cl method [101] and were able 
to reproduce those results; despite a 3% discrepancy in the location of the first resonance, 
near 25 eV, the two methods are in close agreement at the 38.7-eV energy of the FLASH 
measurements. Cl calculations were also performed for the individual dissociation channels, 
and for an isotropic molecular orientation, allowing direct comparison with the FLASH 
cross-section measurement [99]. The Cl calculations give a ratio of ~ 1.7 between the two 
dissociation channels at energies above 35 eV [101], so that the total photoabsorption cross 
section of ~ 0.8 x 10~ 18 cm 2 at 38.7 eV, predicted in the geminal model [67], implies a 
partial cross section of ~ 0.3 x 10 -18 cm 2 to the 4 He + H + channel. The cross-section results, 
shown in Table VII, demonstrate consistency between experiment and theory, although no 
theoretical uncertainties have been assigned and the experimental uncertainty is large. 
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TABLE VII. Photodissociation cross section for 4 HeH + + 7 — > 4 He + H + , from geminal and Cl 
theories as well as from an experiment at FLASH. The geminal result, originally computed for 
both dissociation channels, is corrected for this channel by a factor of 1.7, given by Cl calculations. 



Molecular 

Cross-section 


Orientation 

(10~ 18 cm 2 ) 

Geminal [67] (with Cl [101]) 

Parallel 

~ 0.3 

Cl (adiabatic limit) [101] 

Parallel 

~ 0.46 

FLASH [99] 

Parallel 

0.4(2) 

Cl (adiabatic limit) [101] 

Isotropic 

1.4 

FLASH [99] 

Isotropic 

1.45(7) 


For each event in the FLASH data, the neutral-fragment momentum can be used to re¬ 
construct the initial molecular orientation, under the assumption of fast fragmentation. In 
general, £ — £ transitions peak for orientations parallel to the held, while E — n transitions 
peak when the molecule is oriented perpendicular to the held. For vibrationally cold molec¬ 
ular ions dissociating through the 4 He + H + channel, the measured value of ~ 1 : 3 for the 
E : n contribution ratio [100] agrees reasonably well with the Cl prediction of rv \ ‘ 2 [101], 
There is a clear disagreement in the other channel, however: an experimental measurement 
of E : n ~ 1 : 1, compared to a Cl prediction of ~ 1 : 6. 

Another discrepancy arises in the relative probability of photodissociation to the two 
channels. For vibrationally cold molecular ions, a ratio of <7 He + +H /crHe+H+ = 1-70(48) was 
observed in the later FLASH measurement [100], in agreement with the prediction of about 
1.7 from the Cl method [101]. However, this ratio was found to drop to 0.96(11) when the 
ion beam was not cooled, contradicting the expectation from the Cl potential curves that 
the ratio would rise. 

Without an error estimation from the theory, the significance of these discrepancies be¬ 
tween the Cl model and experiment cannot be evaluated. If the discrepancies hold, they may 
signal the importance of non-adiabatic effects, which were not included in the calculation of 
the Cl potential curves [101]. Such effects are expected to be important to the application 
of the Cl method to the molecular final-state distribution following beta decay in T 2 . 


B. Studies of 3 HeT + and 3 HeH + after beta decay 

1. Instantaneous final-state distribution after beta decay 

In principle, spectroscopy of T 2 gas can be used to measure the instantaneous population 
of accessible 3 HeT + final states after T 2 beta decay, provided that primary radiative tran¬ 
sitions from states excited in beta decay are distinguished from secondary transitions from 
states excited collisionally. One expects that electronic excitations of 3 HeT + will dissociate 
on a time scale of about 10~ 15 s, so any observable radiative transitions must arise from 
excited dissociation products. Consideration of the dissociation channels for each electronic 
excited 3 HeT + state led Jonsell et al. to conclude that only states representing about 16% 
of the total transition probability can result in electronic excited dissociation products that 
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decay via photon emission [54], A calculation of the full probability distribution of dissocia¬ 
tion channels and excitation states is complicated by interference between molecular states 
and has not been attempted. Experimental data on these transitions are sparse: only one 
primary transition has been observed in T 2 spectroscopy, a 468.6-nm line corresponding to 
the 4s —y 3 p transition in 3 He + [102, 103]. 

As seen in Sec. VIA 1, radiative transitions also occur between rotational and vibrational 
levels of 3 HeT + . An infrared emission line (4.69(3) pm) has been observed in T 2 gas and 
identified as the transition between the v — 1 and v = 0 vibrational levels of the 3 HeT + 
electronic ground state [104]. The population of excited rotational and vibrational states 
after T 2 beta decay depends on the beta momentum, but this experiment did not detect 
the beta electrons and was therefore insensitive to this variation. The measured excitation 
probability of the v — 1 level (0.4(2) [104]) thus cannot be compared directly to predictions 
made near the beta endpoint [54], 


2. Branching ratios to electronic excited final states 

The theory can also be probed by measurements of the branching ratios to various regions 
of the final-state spectrum following beta decay in T 2 . A precise measurement of the electron 
energy spectrum about 25 eV below the endpoint would give the branching ratio to the 
electronic excited states of 3 HeT + , which cause a kink in the tritium beta decay spectrum. 
With good energy resolution and a large enough sample window, the change in slope can 
be measured. The energy resolution must be better than 10 eV to resolve the kink, and the 
spectrum must be extended to still lower energies to accurately measure the initial slope. 
Lower energies correspond to much higher rates, imposing a significant additional burden 
on the detector system, and corrections for scattering introduce systematic uncertainty. 

Theory predicts that this branching ratio should be about 43% near the endpoint [54] , but 
no measurement of the branching ratio to electronic excited states has been reported. The 
KATRIN experiment will be able to measure the spectrum in the relevant regime, providing 
the first direct test of the branching ratio to electronic excited states. 


3. First and second moments of FSD from beta decay 

It was pointed out by Staggs et al. [105] that one of the most direct measures of the 
accuracy of the FSD is the comparison of the extrapolated endpoint from beta decay with 
the value expected from mass-spectrometric determinations of the T- 3 He atomic mass dif¬ 
ference, Qa■ If the extrapolated endpoint is obtained from the beta spectrum well below the 
endpoint, it is the average of the individual quantities A k j and differs from the ground-state 
value A 0 o by the first moment of the FSD. Neglecting neutrino mass and the Heaviside func¬ 
tion, which affect the spectrum only at the endpoint, the beta spectrum of Eq. 13 summed 
over final states k becomes 



( 33 ) 
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TABLE VIII. Comparison of zeroth, first, and second moments of theoretical final-state distribu¬ 
tions [10]. 


Reference 

Energy range 

eV 

Efc \ w ko\ 2 

(b(f)k) 

eV 

eV 2 

Fackler et al. [59] 

0 to 165 

0.9949 

-17.71 

611.04 

Saenz et al. [8] 

0 to 240 

0.9988 

-18.41 

694.50 


The summation may be written in terms of binding energies and the atomic mass difference, 

-1 2 


E I’Utt 


(Qa ~b 0 + 2 m e ) 1 - 


Qa ~ b 0 
2M n 


-m e + b [f)k - E, 


E l^wl 2 (<* + &</>* - E ‘Y 


(34) 

(35) 


where terms of order b(f) k m e /M 0 have been dropped and a parameter 6 (the extrapolated 
endpoint energy for zero final-state binding) has been defined for brevity. The summation 
may then be carried out, 


dN _ nv/ry V \PeE e 

M"CF(Z,E e )^- 


1 - ^) o +m 


1 + 


crt 


(5 + (b {f)k ) - E e y 


(36) 


The mean binding energy (b(f) k ) acts as a shift in the extrapolated endpoint S, and the 
variance of = (5 2 ^ fc ) — ( b ^ k ) 2 of the (full) binding-energy distribution enters the expression 
as a shape distortion near the endpoint. Hence, both the first and second moments of the 
final-state distribution can be extracted from data for comparison with theory. Table VIII 
lists the first three moments of the binding-energy distributions for two theories. 

In practice, experiments are not analyzed in this way. Rather, the FSD from theory is 
used to generate the spectrum to be fitted to data, from which values for Qa and rn v can 
be extracted. In addition, only three experiments have used gaseous tritium, and the most 
modern of these (Troitsk [23]) has a scattering contribution to the spectrum at energies 
more than 10 eV below the endpoint. However, the two remaining experiments, LANL [15] 
and LLNL [16] used differential spectrometers and magnetic field configurations designed 
for a broad spectral reach. The two experiments were in good agreement with each other, 
but, as is well known, both found an unexpected excess of events in the endpoint region, 
which is expressed numerically as a negative m 2 . They also yielded concordant values for 
Qa, but only recently has an accurate determination of Qa by a non-beta-decay method, ion 
cyclotron resonance in the Smiletrap apparatus [53] , become available for comparison. Table 
IX shows the results of the LANL and LLNL experiments as originally reported, both having 
been analyzed with the theory of Fackler et al. [59] . The data for those experiments are no 
longer available, but it is possible to estimate the changes that would be produced with the 
use of a more modern theory such as that of Saenz et al. [8] by applying Eqs. 36 and 37. 
The results are shown in the lower half of the table. There is excellent agreement between 
the atomic mass from beta decay and from ion cyclotron resonance, 18589.8(12) eV, and the 
large negative value of m 2 is eliminated in both experiments, subject to the limitations of the 
approximations used. These results provide a striking measure of experimental confirmation 
of the calculations of Saenz et al, especially in the difficult regime of electronic excited 
states. 
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TABLE IX. Atomic mass difference and neutrino mass squared extracted from two experiments, 
in one case with the original 1985 theoretical calculations of the FSD and in the second case with 
a more modern calculation. 



LANL [15] 

LLNL [16] 


As published. Theory: Fackler et al. [59] 


Aoo 

18570.5(20) 

18568.5(20) 

eV 

Qa 

18588.6(20) 

18586.6(25) 

eV 

0 

-147(79) 

-130(25) 

eV 2 

Re-evaluated. Theory: Saenz et al. [8] 


Aoo 

18571.2(20) 

18569.2(20) 

eV 

Qa 

18589.3(20) 

18587.3(25) 

eV 

2 

mi 

20(79) 

37(25) 

eV 2 


4- Branching ratios to molecular and atomic species 

The branching ratio to the bound molecular ion can be extracted from the theory in a 
straightforward way with certain assumptions. Two 1950s mass-spectrometry experiments 
measured this branching ratio for HT [78, 79]; one of the experiments also measured the 
branching ratio for T 2 [79]. The experimental results are consistent with each other but 
disagree starkly with the theoretical prediction. 

Calculations of the dissociation likelihood rely on the theoretical dissociation energy of 
1.897 eV and assume that all electronic excited states are dissociative, i.e. there are no 
fast radiative transitions between the excited states and bound states [54], Under these 
assumptions, and working near the beta endpoint, Jonscll et al. [54] have calculated a 
branching ratio to the bound 3 HeT + molecular ion of 0.39 — 0.57, depending on whether the 
quasibound states above the binding energy dissociate. An absolute uncertainty of 0.2%, 
derived from requiring that the FSD integrate to 100%, is given for calculation of the entire 
spectrum but no explicit uncertainties are indicated for the branching ratios. 

A calculation of the differential spectrum as a function of electron energy would permit a 
more stringent test of the theory than the energy-averaged branching ratio. Experimentally, 
the ability to distinguish between dissociation products ( e.g. between 3 He + +T and 3 He+T + ) 
allows a stronger test than a simple measurement of the dissociation likelihood, yielding 
information about how the electronic states are populated after beta decay. 

The first experimental measurement of molecular dissociation following tritium decay 
was reported for HT by Snell, Pleasanton, and Leming in 1957 [78]. The experiment used 
a mass spectrometer with a conical assembly of ring electrodes that focused ions from an 
equilibrated mixture of HT, T 2 , and H 2 gas into a magnetic analyzer followed by an electron 
multiplier [106]. The measured intensity of the mass-2 peak (H^“) was used to correct the 
other peaks for ionization of the T 2 or HT gas caused by collisions with beta electrons. The 
mass-3 peak (T + or 3 He + ) was corrected for the presence of T 2 in the sample gas, based on 
the ratio of the mass-6 and mass-4 peaks. The correction assumes that HT and T 2 have 
identical dissociation probabilities, which theory does not exclude [54], The final published 
result was a 93.2(19)% branching ratio for HT decay to the bound 3 HeH + ion [78]. 

The following year, Wexler used a mass spectrometer with significantly different ion op¬ 
tics to measure the dissociation probability for both HT and for T 2 [79]. In this apparatus, 






the entire source volume was contained within a cone of ring electrodes, which was followed 
by two distinct deflection stages, one to exclude neutral molecules and one for analysis. A 
measurement with T 2 gas, after correction for an 11.5% HT impurity, yielded a 94.5(6)% 
probability of decay to the bound 3 HeT + . With a pure sample of HT (0.4% T 2 contami¬ 
nation), the probability of decay to the bound 3 HeH + ion was measured at 89.5(11)%, in 
broad agreement (1.2cr) with the Snell et al. measurement [78]. 

In the T 2 dataset, the Wexler apparatus was unable to resolve the difference between 
3 He + + T and 3 He + T + . For an HT source, however, both Wexler [79] and Snell et al. [78] 
found that dissociation into a final state of 3 He + + H was about three times more likely than 
dissociation into 3 He + H + . This is qualitatively similar to the prediction shown in Table VI, 
which yields a ratio of 2.1 for the five electronic excited states considered. 


TABLE X. Branching ratio to the bound molecular ion for HT and T 2 . 


Molecule 

Theory 

Snell et al. 

Wexler 

(Ref. [54]) 

(Ref. [78]) 

(Ref. [79]) 

HT 

0.55-0.57 

0.932(19) 

0.895(11) 

t 2 

0.39-0.57 

- 

0.945(6) 


Table X summarizes theoretical and experimental results for the branching ratio to the 
bound molecular ion. The experimental results for HT and T 2 are in stark disagreement 
with the theoretical predictions. While a problem of this magnitude with the theory seems 
unlikely, it is true that geminal calculations of the bound and continuum states are not done 
in the same basis, and the normalization between the calculations can bias the branching 
ratio. 

To reconcile theory and experiment, other explanations have been advanced for the dis¬ 
crepancy. The applicability of the theory can be questioned in that the experiments inte¬ 
grated over the entire beta spectrum whereas the sudden approximation is valid when the 
electron energy is much larger than atomic binding energies. Another possible mismatch 
between theory and experiment arises from the evolution of the final state before the ions are 
detected. If fast radiative transitions from the electronic excited states to the ground state 
exist, the experimental measurements would have been too slow to prevent repopulation of 
the ground state. At the same time, the measurements may have been too fast for some 
quasi-bound states in the ground-state manifold to dissociate. The time scales for radiative 
decays are, however, expected to be orders of magnitude longer than those for dissociation 
of all but the quasibound states. 

A number of experimental issues have also been identified. The experiments may not have 
properly accounted for contamination of the mass-6 signal by produced via ionization, 
artificially inflating the measured branching ratio to the bound molecular ion. This risk was 
not unknown to the experimenters, who took steps to mitigate it. 

Wexler himself favors the explanation that the relative efficiencies between ion species 
were poorly understood, as the acceptance of both mass spectrometers depended strongly 
on the initial transverse energy of the ion [7, 54, 79]. This transverse energy is dependent 
on the ion species and can range up to tens of eV following dissociation of excited states of 
3 HeT + , although most of the dissociation processes should lead to ions in the energy range 3 
- 13 eV. As computed in Sec. V C, the ion energies resulting from excited-state dissociation 
are larger than the ~ 1-eV energies for mass-3 fragments in the breakup of the ground state, 
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but whether this accounts for the experimental results is not possible to determine without a 
model for the acceptance of the mass spectrometers. A more telling observation, however, is 
that in the decay of HT the energies of the mass-3 fragments are lower than in the decay of 
T 2 . That is consistent with Wexlcr’s suggestion because the measured branch to the bound 
final state HeH + is smaller than that to HeT + , perhaps due to better efficiency for detecting 
the dissociation fragments. One may also surmise that while dissociation is energetically 
allowed from the ground-state manifold above 1.897 eV excitation, it is strongly hindered by 
the angular momentum barrier. A much larger fraction of the HeT + ground-state manifold 
can potentially decay this way than for HeH + , and yet the data show the opposite behavior. 

The disagreement between theory and experiment has not been satisfactorily explained, 
although many sources of possible unquantified experimental error have been proposed. No 
data are available to test these explanations, however. Further measurements with the 
potential to resolve this tension are desirable. 


C. Desiderata for a modern experiment 

A modern dissociation experiment could more closely reproduce the conditions for which 
the calculations are performed. Detecting the ion in coincidence with a beta electron of 
measured energy would allow the experimenter to examine the specific regime where the 
sudden approximation is valid and to study the variation of the dissociation fraction with 
electron energy. The acceptance of the instrument for ions with a range of initial kinetic 
energies needs to be quantifiable. Measurement of the ion energy distribution would provide 
a stronger test of the model. Complementary information is also available in the coincident 
photon spectrum but the expected emission falls in the vacuum ultraviolet regime, making it 
difficult to instrument. Operating conditions must be such that charge exchange is a minor 
and quantifiable perturbation. 

A way of implementing many of these objectives is the use of semiconductor detectors 
and low-pressure tritium in uniform, coaxial electrostatic and magnetic fields. Mass separa¬ 
tion is achieved by time of flight, and the held arrangement offers high efficiency. When the 
magnetic held strength is sufficient to collect ions regardless of their transverse momentum, 
the species-dependent efficiency changes can be eliminated. The radial excursions of the 
ions can, moreover, be mapped to provide information about their energies and to provide 
assurance that all have been detected. Higher detection efficiency allows the source pressure 
to be lowered, reducing charge exchange, which can artificially lower the measured dissoci¬ 
ation probability. An experiment utilizing this approach could more closely reproduce the 
conditions of the calculations and provide a direct test of specihc aspects relevant to the neu¬ 
trino mass measurement. Such an experiment, the Tritium Recoil-Ion Mass Spectrometer 
(TRIMS), is under construction at the University of Washington. 


VII. DISCUSSION AND CONCLUSIONS 

A. Impact on tritium neutrino mass experiments 

In this section we aggregate and, where possible, quantify the various ways in which FSD 
uncertainties contribute when a gaseous tritium source is used to measure neutrino mass. 
These fall into 3 groups: theoretical uncertainties in the FSD itself, uncertainties in the 
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degree of temperature equilibration for T 2 , and uncertainties in the isotopic composition of 
the source gas. 

The KATRIN experiment has sufficient statistical power that data-taking can be con¬ 
centrated in the last 20 eV of the spectrum, which removes the theoretical uncertainties in 
electronic excitation of the molecule as a major concern. There is remaining uncertainty 
in the width of the ground-state manifold of rotational and vibrational excitations, but we 
have shown that the broadening has a very simple origin, mainly zero-point motion. Indeed, 
the semiclassically derived analytic expression yields a variance that agrees with the full 
theoretical calculation to 7%. Beyond this, a quantitative uncertainty estimate is lacking, 
and knowledge of the variance at the 1% level has been assumed in the design of experiments 
like KATRIN. We have reviewed a variety of tests of the theory, finding generally excellent 
agreement, with the one serious exception being the branching ratio to the bound mass-6 
ground state manifold. A new experiment would provide substance for a re-evaluation of 
the theoretical uncertainties. 

An accurate characterization of the composition of the source is necessary for KATRIN. 
The source gas is high-purity T 2 . To determine the isotopic composition, the KATRIN 
collaboration has developed a laser Raman spectroscopy system called LARA. This system 
has achieved 0.1% precision [38] and better than 10% accuracy [39] in measurements of 
the isotopic composition. In principle, a laser Raman system can also provide information 
about the ortho-para ratio. However, due to the difficulty of in situ measurement, LARA is 
located at a high-pressure stage prior to cooling and injection into the source. The KATRIN 
collaboration is studying an extension of the LARA system to measure the ortho-para ratio 
and is conducting ongoing simulation work on the evolution of the ortho-para ratio and 
other source parameters. 

The KATRIN wiudowless, gaseous tritium source vessel will be maintained at a temper¬ 
ature of 30 K. In thermal equilibrium at this temperature the ortho-para ratio is approxi¬ 
mately 1:1 and states with J > 1 are not appreciably populated. The time each molecule 
spends in the cooled source, however, is short compared to the spin relaxation time. The 
ortho-para ratio of the gas within the source is therefore expected to be close to 3:1. 

Disequilibrium in the source is not confined to the ortho-para ratio because depopulation 
of higher excited states in free space requires quadrupole transitions that are very slow. 
The de-excitation process is therefore predominantly collisional and apparatus-dependent. 
Incomplete thermalization of these excited states would be a source of uncertainty if undi¬ 
agnosed. 

These sources of uncertainty in the FSD translate directly to an uncertainty in the neu¬ 
trino mass-squared. Robertson and Knapp [10] have shown that any neglected contribution 
to the variance of the FSD, Acr| SD , modifies the extracted neutrino mass-squared by 

Ami - —2AcTp SD . (37) 

Doss et al. [9] calculated the final state distributions arising from the lowest four rota¬ 
tional states of T 2 and the lowest two states of HT and DT, i.e. those populated in a 30-K 
thermal source. The FSDs were binned with 0.01-eV resolution compared to the 0.1-eV 
resolution used in reporting previous results [8]. We have estimated the variance of each 
binned distribution in two ways: using the central bin energy value and the reported mean 
energy value. We take the average of the two results as the best estimate of the variance and 
half the difference as the width (standard deviation) of the error distribution. The mean ex¬ 
citation energies and estimated variances of the FSDs are listed in Table XI. Unfortunately 
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TABLE XI. Mean excitation energy and variances extracted from the FSD calculations of refer¬ 
ence [49]. There is a small contribution to the variance (< 0.004 eV 2 ) from binning. 


Source 

J 

Mean E exc (eV) 

<6 (eV 2 ) 


0 

1.752 

0.194 

t 2 

1 

1.751 

0.206 

2 

1.750 

0.215 


3 

1.749 

0.262 

DT 

0 

1.752* 

0.175 

1 

1.752* 

0.188 


*Shifted to compensate for different recoil kinetic energy [49]. 


the distributions for higher rotational states of T 2 were not available, and distributions for 
HT are not available with the required binning resolution. Future calculations of the FSD, 
such as calculations using the configuration-interaction method, will be helpful in expanding 
and improving the estimates of the variances. 

Figure 7 compares the semi classical variances calculated for initial states (0, J ) in T 2 
using Eq. 30 to the variances estimated from the calculations of Doss et al. [ 9 ]. From 
the figure we conclude that the semiclassical model is a good proxy for the FSD variance. 
The difference between the two is about 7% and independent of J. Of this difference, 1% 
is attributable to our more accurate result for because all contributions to the 

I cCjillcLX. 

variance are proportional to p 2 /2m. Given the limited set of full FSD calculations available, 
we use the semiclassical variances to estimate the systematic errors associated with various 
experimental parameters. 

After shifting the excitation energy to compensate for differences in the recoil kinetic 
energy, the effective mean excitation energy of each of the FSDs corresponds to the same 
laboratory endpoint energy for each isotopolog. Thus the variance of the summed distribu¬ 
tion can be taken as the sum of the variances for each isotopolog i and each rotational state 
J, weighted according to their populations fi and Pj ti for isotope and rotational state, respec¬ 
tively. An additional variance contribution arises from the translational Doppler broadening 
<7 2 rans at a given temperature T. The overall variance a 2 ot of the line broadening can be 
derived: 

2 = pf_ \2fH p 2a 2 m 2 J(J + 1) 
aj ’ 1 2m 3 m zp ^ ?>R$m 

°FSD,i = 

J 

2 p 2 2mksT 

<Arans,i = ^ ^ + m 

a tot = fi (^FSD,i + ^trans,*) 

i 

The Pjj weights are given by a Boltzmann distribution for the temperature T. (The trans¬ 
lational and rotational temperatures need not be the same). The probability distribution 
is calculated independently for each isotopolog and summed according to the activity frac¬ 
tion fi of each isotopolog in the source. The source activity may be expressed in terms of 


(38) 

(39) 

(40) 

(41) 
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FIG. 7. Comparison of the variance of the ground-state-manifold FSD produced in T 2 decay as 
calculated in the semiclassical model, Eq. 30 (solid curve, red online), with variances taken from 
calculations for states up to J = 3 described in Ref. [9] (blue dots). 


a parameter ex that is the equivalent fractional activity of the gas compared to pure T 2 . 
Additionally the ratio of HT to DT in the source gas n = /ht//dt is used to characterize 
the makeup of the active contaminants. Eq. 42 shows the functional form of the isotopic 
weights. 

( 2ex — 1 ,i — T 2 

fi= < 2(l-e T )/(l + «) ,i = DT (42) 

[ 2(1 — 6t)k/(1 + a) ,i = HT 

Neglecting inert isotopologs H 2 , HD, and D 2 , ex is confined to the range 0.5 < ex < 1 and is 
assigned a reference value of 0.95 as in the KATRIN Design Report [6]. The reference value 
of k is taken to be 0.1 because the fractional distillation process results in higher levels of 
deuterium than of protium. 

Table XII shows the rotational-state distributions for T 2 thermal 30 K, thermal 300 K, 
and nonthermal 30 K (A = 0.75) sources along with the semiclassical FSD variances. Also 
shown is the contribution each state makes to the total FSD variance of the source in each 
configuration. The rotational-state distributions for DT and HT are shown in Tables XIII 
and XIV, respectively. (The rotational-state energies differ slightly from those given by 
Doss [9], possibly because centrifugal stretching is not included here.) The rotational states 
up to J = 7 contribute significantly at room temperature and further work is necessary 
to provide an accurate assessment of the systematic error associated with the experimental 
uncertainty in the rotational-state distribution. Measurement of the rotational-state tem¬ 
perature and calculations of the higher rotational-state FSDs would significantly improve 
the error estimates. 

To quantify the impact of using an incorrect FSD to analyze neutrino-mass data we ex- 


33 



TABLE XII. Rotational-state distributions for To at 30 K and 300 K. The energies are those used 
in Ref. [49] and variances are from the semiclassical width, Eq. 30. Probabilities P are calculated 
from the partition function (Eq. 16) using the energies listed in the table and the contributions to 
the total FSD variance are computed accordingly. 


J 

Ej 

(meV) 

a J, t 2 

0 

CO 

Thermal 

300 K, 

Thermal 

30 K, 

A = 0.75 

(eV 2 ) 

P (%) Var Contr 

P (%) Var Contr P (%) Var Contr 

0 

0.00 

0.1762 

43.70 

0.0768 

4.73 

0.0083 

24.6 

0.0434 

1 

5.01 

0.1875 

55.70 

0.1040 

35.00 

0.0656 

75.0 

0.1410 

2 

15.02 

0.2103 

0.62 

0.0013 

13.20 

0.0277 

0.35 

0.0007 

3 

30.05 

0.2445 

0.01 

0.0000 

30.70 

0.0752 

0.01 

0.0000 

4 

50.08 

0.2900 

0.00 

0.0000 

6.03 

0.0175 

0.00 

0.0000 

5 

75.11 

0.3469 

0.00 

0.0000 

8.33 

0.0289 

0.00 

0.0000 

6 

105.16 

0.4152 

0.00 

0.0000 

1.02 

0.0042 

0.00 

0.0000 

7 

140.21 

0.4949 

0.00 

0.0000 

0.90 

0.0045 

0.00 

0.0000 

8 

180.27 

0.5859 

0.00 

0.0000 

0.07 

0.0004 

0.00 

0.0000 

9 

225.34 

0.6883 

0.00 

0.0000 

0.04 

0.0003 

0.00 

0.0000 

10 

275.42 

0.8022 

0.00 

0.0000 

0.00 

0.0000 

0.00 

0.0000 



FSD Variance 

0.1830 


0.2330 


0.1850 


amine the differences in variances that arise due to changes in temperature, isotopic purity 
and ortho-para conditions. For small deviations from the operating parameters the corre¬ 
sponding error in the extracted neutrino mass-squared can be derived from Eq. 37. Below, 
we derive the functional form for ortho-para ratio errors, temperature fluctuations and errors 
in the isotopic composition. The nominal source parameters are shown in Table XV. 

The temperature of the source is a key parameter determining the width of the final state 
distribution. As previously stated, the rotational states of homonuclear T 2 do not equilibrate 
on short time scales [107] and the exact time required for thermalization in the KATRIN 
source depends not only on the gas density but also on the materials the gas contacts (i.e. 
walls, permeators, etc.). The temperature changes the initial rotational-state distribution of 
the source as seen in the partition function. For small fractional changes in temperature the 
exponential factors can be expanded, and the resulting shift in variance can be expressed 
in terms of the fractional temperature change. For a cryogenic source only the J — 0 and 
J = 1 states contribute significantly and the shift in FSD variance for a given isotopolog 
simplifies to a single term, which may be written: 


^FSD.i 




Ej ti - E nji ST 

T 


kT 


P^Pi-P- 2 (ol-al)5T. 


(43) 

(44) 


Table XVI shows the translational Doppler variance temperature-variation coefficients for 
T 2 , DT and HT, computed from: 


da, 


p 2 2 mk 


B 


trans.z 


2 m m s i + m 


ST. 


(45) 
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TABLE XIII. Rotational-state distributions for DT at 30 K and 300 K. The energies and variances 
are from the semiclassical model (see Eq. 30). Probabilities are calculated from the partition 
function (Eq. 16) using the energies listed in the table and the contributions to the total FSD 
variance are computed accordingly. 


J 

Ej 

(meV) 

2 

a J, DT 

(eV 2 ) 

30 K, 

Thermal 

300 K, 

Thermal 

P (%) Var Contr 

P (%) Var Contr 

0 

0.00 

0.1578 

78.70 

0.1242 

11.61 

0.0183 

1 

6.25 

0.1692 

21.02 

0.0356 

27.36 

0.0463 

2 

18.76 

0.1919 

0.28 

0.0005 

28.11 

0.0540 

3 

37.52 

0.2261 

0.00 

0.0000 

19.05 

0.0431 

4 

62.53 

0.2716 

0.00 

0.0000 

9.31 

0.0253 

5 

93.80 

0.3285 

0.00 

0.0000 

3.39 

0.0111 

6 

131.32 

0.3968 

0.00 

0.0000 

0.94 

0.0037 

7 

175.09 

0.4765 

0.00 

0.0000 

0.20 

0.0010 

8 

225.12 

0.5675 

0.00 

0.0000 

0.03 

0.0002 

9 

281.40 

0.6700 

0.00 

0.0000 

0.00 

0.0000 

10 

343.93 

0.7838 

0.00 

0.0000 

0.00 

0.0000 



FSD Variance 

0.1603 


0.2029 


The shifts in variance due to the FSD and translational effects are additive, and each 
contributes to the overall shift in the extracted neutrino mass-squared according to Eq. 37. A 
temperature change of 0.15 K from the nominal 30 K results in a shift in extracted neutrino 
mass-squared of 0.11 x 10~ 3 eV 2 . 

In reality both thermal fluctuations and inaccuracy in the measurement of the temper¬ 
ature contribute to the uncertainty on neutrino mass. It is reasonable to assume these 
are uncorrelated errors and thus two independent thermal factors appear in the error bud¬ 
get. The expected temperature fluctuations and uncertainties are taken from the work of 
Grohmann et al. [37, 108]. 

The isotopic purity of the source plays a major role in neutrino-mass experiments because 
the width of the FSD varies significantly between isotopologs. In addition to the dependence 
on the tritium activity fraction ex, there is a dependence on the relative population k of 
contaminants HT and DT. Tables XII, XIII and XIV show the variance of the distribution 
for 30-K sources of tritium-containing isotopologs. The T 2 results include the thermal source 
as well as the nonthermal source with A = 0.75. The large differences in the FSD variances 
between HT, DT and T 2 demonstrate the importance of knowing the isotopic purity. The 
shift in the variance that occurs when the tritium purity of the source 6 t changes can be 
written 


(5a 2 



2 

1 + K 


a 


DT 


2 K 

1 + K 



fit T, 


(46) 


where of is the sum of the FSD (Eq. 39) and translational (40) terms. Similarly, the 
dependence on k takes the form 


Sa 2 


2(1 — e T ) 
(1 + k ) 2 


a 


DT 


+ a 


HT 


Sk. 


(47) 
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TABLE XIV. Rotational-state distributions for HT at 30 K and 300 K. The energies and variances 
are from the semiclassical model (see Eq. 30). Probabilities are calculated from the partition 
function (Eq. 16) using the energies listed in the table and the contributions to the total FSD 
variance are computed accordingly. 


J 

Ej 

(meV) 

2 

a J,WT 

(eV 2 ) 

30 K, 

Thermal 

300 K, Thermal 

P (%) Var Contr 

P (%) Var Contr 

HT 0 

0.00 

0.1251 

94.09 

0.1177 

18.12 

0.0227 

1 

10.00 

0.1365 

5.91 

0.0081 

36.93 

0.0504 

2 

29.99 

0.1592 

0.00 

0.0000 

28.40 

0.0452 

3 

59.98 

0.1934 

0.00 

0.0000 

12.46 

0.0241 

4 

99.97 

0.2389 

0.00 

0.0000 

3.41 

0.0082 

5 

149.95 

0.2958 

0.00 

0.0000 

0.60 

0.0018 

6 

209.94 

0.3641 

0.00 

0.0000 

0.07 

0.0003 

7 

279.91 

0.4438 

0.00 

0.0000 

0.01 

0.0000 

8 

359.89 

0.5348 

0.00 

0.0000 

0.00 

0.0000 

9 

449.86 

0.6373 

0.00 

0.0000 

0.00 

0.0000 

10 

549.83 

0.7511 

0.00 

0.0000 

0.00 

0.0000 



FSD Variance 

0.1258 


0.1526 


TABLE XV. Reference values of parameters used in estimating FSD and Doppler contributions to 
the projected uncertainty in the extracted m 2 for KATRIN. 


Parameter 

Value 

Source temperature 

T = 

30 K 

Ortho fraction 

A = 

0.75 

Tritium fraction in WGTS 

Ct = 

= 0.95 

Ratio of DT to HT 

K = 

0.1 


Starting from the nominal source parameters (Table XV) and introducing an uncertainty 
of 1% on the atomic purity of the source would lead to a uncertainty on the neutrino mass- 
squared of 0.96 x 1CT 3 eV 2 . While conflicting previous results have led to confusion over 
the impact of errors in the measurement of isotopic purity [9, 39], our results are consistent 
with the earlier published work of Doss et al. [9] which concluded that it plays a major role. 

The impact of the ortho-para condition of the source can also be derived from Eq. 41 
by considering a slight reordering of rotational states. Due to the two-state nature of the 


TABLE XVI. Variation with temperature of the translational Doppler contribution to the variance 
for a source near 30 K, calculated from Eq. 45. 


Source 

Str jffT (10 -3 eV 2 /K) 

t 2 

0.147 

DT 

0.176 

HT 

0.220 
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homonuclear system, the state distribution for T 2 is often separated out in terms of the even 
(para) and odd (ortho) states. The sum of probabilites for all the odd states is the ortho 
fraction of the source: 

A = Yi Pj. (48) 

J odd 

The variances of the ortho and para states can then be considered separately and even 
normalized independently to yield ortho and para state probabilities, labeled P 0 rtho,j and 
Tpara,j respectively. The total variance is then the sum of two states weighted according to 
the A factor. 


^FSD^ — A ^ ^ -Portho,J cr J + (1 A ) ^ ^ -Ppara.J^J (49) 

J odd J even 

= Acr ortho + (1 _ '^) Cr para- (50) 

If the probabilities within the ortho (para) state relative to the other states are not 
changing then the impact of the ortho-para transitions can be assessed in terms of the 
independent ortho and para state variances. Under this assumption, the dependence on 5 A 
is simply characterized by the difference in the FSD variances arising from the ortho and 
para distributions: 

A(T FSD = (Cktho — °para)^ A - (51) 

For cryogenic sources the equation of the shift in neutrino mass-squared further simplifies, 
only depending on the difference in the variances of the J = 0 and J = 1 states. For small 
changes in temperature which do not appreciably change the occupation of the higher states, 
the shift in variance is independent of temperature. The contributions from DT and HT 
remain unchanged as ortho-para considerations only apply to the homonuclear isotopolog. 
The effect of a change in ortho-para ratio on the extracted neutrino mass-squared is given 
by: 

\Aml\ ~ 2(2e T - l)(oj=i “ (T j=o) (5A - (52) 

Given the relatively short time that molecules will spend at cryogenic temperatures in 
the KATRIN source, the ortho fraction is expected to be close to 0.75, corresponding to 
the 700 K permeator through which the gas passes in atomic form. A lower bound of 0.57 
is set by the beam-tube temperature of 30 K. If A lies at an unknown value between these 
bounds the corresponding uncertainty on the extracted neutrino mass-squared would be 
3.8 x 10 -3 eV 2 . Fortunately this is not expected to be the case and early simulations indicate 
that even in pessimistic scenarios only 3% of the ortho source molecules will transition 
from the ortho state to the para state [109]. These KATRIN simulations show a shift in 
neutrino mass-squared of 0.48(7) x 10 -3 eV 2 due to ortho-para transitions. Our calculation is 
0.44 x 10~ 3 eV 2 , in good agreement with the results of the simulation. Thus under standard 
scenarios the ortho-para ratio is not expected to contribute significantly to the uncertainty 
on the neutrino mass-squared. 

While not considered a significant concern for KATRIN, from an experimental perspective 
the ortho-para ratio warrants more study as the A factor and associated systematic error 
can potentially be measured. Current work by the LARA subgroup of KATRIN focuses on 
how to measure the ortho-para ratio using a modified version of the setup used to measure 
the isotopic ratio. 

Table XVII summarizes the projected role of molecular effects on the KATRIN measure¬ 
ment for selected reference values of parameters, showing the sources of systematic error 
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TABLE XVII. Summary of molecular-related sources of systematic shift in extracted neutrino mass- 
squared, the projected accuracy on the experimental parameters and the individual effect on m 2 
for the nominal KATRIN parameters shown in Table XV. The accuracy of theoretical calculations 
of the width is taken as 1% in accordance with the KATRIN Design Report [6] but further study is 
necessary to validate this number as discussed in the text. The achievable experimental uncertainty 
on the rotational-state temperature is being studied but is not known at this time. 


Source of systematic shift 

Target accuracy 

O‘syst(m^)[10 3 eV 2 ] 

FSD theoretical calculations 

Aofsd/cfsd —1% 

6 

temperature calibration 

AT /T\ < 0.005 


- translational 


0.05 

- FSD 


0.06 

temperature fluctuations 

AT/T < 0.001 


- translational 


0.009 

- FSD 


0.01 

ortho-para ratio 

AA/A < 0.03 

0.44 

isotopic impurities 
- tritium purity 

|Act/ctI A 0.03 

2.9 

- ratio of HT to DT 

Ak/k = 0.1 

0.03 

higher rotational states 

^ ^rotational/^rotational | 

= 0.1 1 


associated with molecular excitations, the projected accuracy on the parameters and the 
corresponding systematic error on the neutrino mass-squared. 


B. Summary 

The use of molecular tritium in experiments to measure the mass of the neutrino neces¬ 
sitates a quantitative understanding of the role of molecular excitations in modifying the 
shape of the observed beta spectrum in the vicinity of the endpoint. Electronic excita¬ 
tions are important but as experimental sensitivity has improved, the focus has increasingly 
shifted to the rotational and vibrational excitations of the daughter molecule in its elec¬ 
tronic ground state. Those excitations modify the spectrum at the endpoint, whereas the 
electronic excitations set in some 20 eV below the endpoint. The KATRIN experiment, by 
virtue of its high statistical sensitivity and excellent resolution, will be able to concentrate 
its data-taking in the last 20 eV of the spectrum. 

Detailed quantum calculations of the molecular final-state spectrum have been published, 
and will be used in the analysis of forthcoming experiments. We have shown that the ground- 
state rotational and vibrational manifold is fundamentally a Gaussian distribution with a 
variance determined almost completely by zero-point motion of the nuclei in the parent 
molecule. Structure is imposed on that smooth distribution by the quantized nature of 
the spectrum of hnal states. The simplicity of the underlying mechanism suggests that 
the theoretical prediction of the width of the ground-state manifold should indeed be very 
reliable, as has been assumed in the design of experiments such as KATRIN. Calculations 
using the configuration-interaction method would provide independent uncertainty estimates 
as well as a comparison to the geminal method calculations. This would be a significant 
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improvement over the current assessment of errors, which is based solely on the integral of 
the entire spectrum. 

Thermal excitations of rotational states play a major role for the homonuclear molecule 
T 2 since equilibration of the ortho-para ratio is not immediate. The contribution to the 
width of the ground-state manifold from rotational-state excitations is relatively small if 
the molecule is in thermal equilibrium at a temperature near 30 K, but is significant if the 
distribution remains effectively at 300 K because of the slow thermalization of the ortho¬ 
para systems. Thus the ortho-para ratio must be determined by design or diagnosis. There 
is a need for additional theoretical calculations to map out the contributions of states with 
J > 4. These issues could be circumvented in an experiment that uses HT instead of T 2 . 
Another advantage of using HT is that at 30 K the final-state distribution variance (in 
the ground-state manifold) is 2/3 as large as it is in T 2 . These advantages are somewhat 
counterbalanced by a loss of statistical power caused by the dilution of the activity by 
protium and by the lower source column density caused by the lower molecular mass. 

Although no means is known for a direct experimental measurement of the final-state 
energy spectrum (other than beta decay itself), the theory makes numerous testable predic¬ 
tions. The energies of states in the ground-state manifold are in very precise agreement with 
theory. Re-evaluating the analysis of the Los Alamos and Livermore gaseous tritium exper¬ 
iments with the current theoretical model produces excellent agreement between the atomic 
mass difference determined by beta decay and by ion cyclotron resonance. Furthermore, it 
eliminates the large negative values of rri'l originally reported in those experiments. 

On the other hand, the measured branching ratios to the bound molecular ions 3 HeT + 
and 3 HeH + are in the range 90-95%, in strong disagreement with the theoretical prediction 
of 39-57%. This discrepancy has endured for more than 50 years and a number of possible 
explanations for it have been suggested. Several avenues are now open for progress toward a 
resolution. New work with the configuration-interaction method is underway [68] and may 
result in the first independent theoretical cross-check of modern calculations in the geminal 
basis. A new, direct measurement with beta-ion coincidence information is now feasible with 
modern instrumentation and is being attempted. We have presented schematic calculations 
of the recoil-fragment energy spectra following dissociation, a new and potentially testable 
aspect of the theory. Finally, the KATRIN experiment itself will be able to determine the 
relative fraction of population of the electronic ground and excited states. With a theoretical 
cross-check, new experimental information, and insight into the basic mechanism for final- 
state broadening, one can anticipate increased confidence in quantifying the role it plays 
when extracting a value for the neutrino mass from data. 
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